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ABSTRACT 

An On-Line Modeling System (OLMS) which should be 
particularly useful to the biomedical community is described which 
utilizes a process called "overlaying” to simplify loading and 
executing programs. OLMS has been constructed and implemented as an 
interpreter; that is, a program that translates and then executes 
each source statement in sequence where these two operations follow 
each other in close time proximity. It has been written so that it 
can be run under the IBM S/360 and the PDP-10 operating systems, but 
has been implemented only on the PDP-10. OLMS operates from a command 
language that is executed inter pretively through a set of closed 
subroutines. The user may run a job, save data, get data, and display 
data. While his program is running, he can access the whole memory. 
This report consists of a description of OLMS and programmer's 
reference manuals for both the OLMS graphical system and OLMS. (JYj 
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The Division of Computer Research and Technology, NIH, will 
issue on an irregular basis technical documents which we believe 
will be of particular interest to the biomedical community. 

These reports will include detailed descriptions of relevant 
computer programs and instructions in their use (as well as some 
theoretical background) , in hopes that interested scientists will 
be encouraged to gain first-hand experience in applying them. 

In some cases, such reports may serve as foci around which DCRT 
will structure training courses to expand the knowledge and 
experience of NIH staff in applying computer science to problems 
of research and management. Circulation of these reports within 
the biomedical community broadly is, of course, encouraged. 




A. W. Pratt, M.D., Director, DCRT 
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PREFA CE 

This report is intended to introduce a simple modeling system, that 
is easy to use, flexible and general enough to allow a scientist or 
programmer to extend it using Fortran IV. 

There are three parts to this report; 

Part I : The On Line Modeling System , 

Part II : A Programmers' Reference Manual For The On Line 

Modeling Graphical System, and 

Part III : A Programmers' Reference Manual For The On Line 

Modeling System. 

Before using this system a user should read all of Part I and at least the 
miscellaneous Section in Part III. 

i 
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SECTION I 



I 

1.1 INTRODUCTION 



Access to computers is limited by many constraints; namely, memory 
size, memory cost, computer structure, computer language, internal 
operating systems and monitors, etc. The design and implementation of 
any program that a user wishes to interact with is directly influenced 
by the above factors and perhaps others that are not listed. Most 
computers are designed for general purpose usage and as such the 
generality given any particular user is a compromise between the speed 
of loading and executing programs. Mathematical models of complex 
processes require complicated programs which tend to grow in size as 
the problem solving capability increases, and new knowledge of the 
system being modeled is incorporated. This concept is clear if we 
think of modeling a physiological structure. Independent unit structures 
may be modeled, and at some point, connected together as a model for 
a complete structure. An example is given in Section IV showing how the 
input output relations of the Hodgkin-Huxley Nerve Equations can be 
studied. 

Since most models grow as the investigator proceeds, the total storage 
requirements for instructions in a given program may exceed the available 
storage capacity of the machine. It then becomes necessary to use 
repeatedly the same blocks of the computer's internal storage in order 
to execute the total program. The process whereby data or subprograms 
replace other data or subprograms which are no longer needed, is called 
overlaying . This overlaying process provides a more efficient utilization 
| of memory capacity and increases the flexibility of the internal programming 
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structure. The On-Line Modeling System discussed here utilizes this over- 
lay concept to simplify the process of loading and executing programs. 

A program that translates and then executes each source statement 
(e.g., an individual Fortran statement such as x = y) in sequence, where 
these two operations follow each other in close time proximity is called 
an interpreter . The use of an interpreter is advantageous if rapid 
response to modifications of the source statements is required. For this 
reason The On-Line Modeling System has been constructed and implemented 
as an interpreter. 

Present day operating systems and monitors are influenced by the needs 
of many people; however, they may not completely fulfill the needs of any 
particular user. Furthermore, many generalized operating systems are 
designed in such a way that they cannot normally allow the user to run 
programs, save data, get data, and display data by pushing buttons on a 
keyboard. The On-Line Modeling System described in this report has been 
written so that it can be run under the existing (S/360) and (PDP-10) 

Operating Systems. However, it has been implemented only on the PDP-10. 

The OLMS operates from a command language that is executed interpret ively, 
through a set of closed subroutines . This command language has eight modes 
of operation which are discussed later in this report. In one of these, 
the Function Key Mode, the system can be operated by a push button key- 
board, allowing a user to manipulate the OLMS without having to know the 
details of the command language. 



P In summary, the OLMS is a system that allows a user to run a job, save 

data, get data, generate data, and display data. When the user's program 
is running he can access the whole memory (i.e., the OLMS is not in memory 
when the user's program is running). This is accomplished by using the 
system overlay loader to overlay the OLMS. The interpreter for the OLMS 
command language gives the user flexibility in using the various modes of 
the system. In the Function Key Mode a user can interact with the system 
by pushing buttons and answering questions. As the user answers the 
questions a syntactic command is generated and passed to the interpreter 
to test its acceptance. The advantages of such a system are as follows: 

1) A user can use the system without any display programming 
or any prior knowledge of how the display works. 

2) A user can run a program and see immediate relationships 

[ among the variables (i.e., graphs can be generated by 

pushing buttons) . 

3) When the user's program is running he has the capability 
of using the central processor, fast memory, and all I/O 
devices . 

Finally, the handling of programs that are too large to fit in memory 
is possible since the files are each organized and the user can interact 
with the total organization. Therefore, the OLMS provides investigators 
who must use large programs and large amounts of experimental data with a 
flexible and responsive .omputer facility which allows them to explore all 
the characteristics of their data. Furthermore, the individual user may 
modify the OLMS with some prior knowledge of Fortran. 

ft 
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1.2 Example Problem 



This example shows how the OEMS can be used to study the input and 

2 

output relations of a simple quadratic equation of the form Y = AX + BX + C 
The algor ithm listed below shows how this equation can be implemented on the 
OLMS. 



Algorithm 

The system buffer is defined as a named common bloc’:.; namely, COMMON/ COMBUF/ 
SYSBUF(384) . 

Al. [Input coefficients by teletype] 

A,B,C-< INPUT 

A2. [Input the number of points to generate by teletype]. 

NPTS* INPUT 

if (NPTS>128) then go to A9 
A3. [Initialize I], 

A4. [Generate points]. 

if (I>NPTS) then go to A8. 

A5 . [Generate X value] . 

SYSBUF(I)-< 1 



1-4 



A6 . [Generate Y value] . 

SYSBUF(128+I) = A* (1**2) + B*I + C 

A7 [Increment I] . 

I+---I+1 
Go to A4 

A8 . [Write SYSBUF on disk] . 

CALL WRR(l) 

A9 . [Load Interpreter] . 

CALL INTERP(l) 

A10 . Terminate 



In this algorithm the routine WRR(l) writes the content of the system 
buffer on the disk. The routine INTERP(l) is a call to the Loader to load 
the OIMS interpreter. 

A Fortran IV program is given to show an implementation of this algorithm. 
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FORTRAN IV PROGRAM 



COMMON/ COMBUF/SYSBUF (384) 

TYPE 12 

ACCEPT 1,A,B,C 

1 FORMAT (3F) 

TYPE 13 
ACCEPT 2 ,NPTS 

2 FORMAT (I) 

IF(NPTS.GT. 128) TYPE 11 
IF(NPTS.GT.128)GO TO 9 
I = 1 

4 IF(I.GT.NPTS)GO TO 8 
SYSBUF(I)+I 

SY SBUF (128+1) =A* FLOAT ( I ) * * 2+B*FLOAT ( I ) +C 
I = 1+1 
GO TO 4 

8 CALL WRR(l) 

9 CALL INTERP(l) 

11 FORMAT ( ' ERROR NPTS IS, GREATER THAN 128 7 s ’) 

12 FORMAT ( ' INPUT A,B,C'/' ') 

13 FORMAT ( ' INPUT NPTS'/' ') 

END 



Steps used to get a job ready for the OLMS. 

(1) Place cards on disk 
R PIP 

DSK.-NAME CDR: 

(.2) Compile job 
COM NAME 

(3) Load job 

LOAD NAME, INTOP 

(4) Save job 

SA DSK:NAME 
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The information needed to operate the OLMS can be acquired by reading 
The Function Key Mode in Section III, Section V, and the following 
semantics in Section II of this report: 

4 . 1 (Save Statement) 

5 . 1 (Get Statement) 

6.1 (Delete Statement) 

9.1 (Clear Statement) 

10.1 (Rename Statement) 

12.1 (Plant Statement) 

16.1 (Start Statement) 

18.1 (Kill Statement) 

19.1 (Enter Statement) 

20.1 (Display Statement) 

22.1 (Dout Statement) 

23.1 (Top Statement) 

25.1 (Transfer Statement). 

Now you are ready to use the OLMS; therefore, branch to the computer, other- 
wise, proceed to add depth to your understanding. 



1.3 Generalized Statement of Problem 



Given that a program is too large to fit in the existing computer 
memory, we may divide it into Ch parts (where j = l,2,...,n). The sub- 
division is made so that the user may interact with any . 

Each Cj is loaded for execution with an overlay structure. Suppose 
that the part has been loaded and executed. Depending on the output 
from the user may want to execute again or one of the other 
C 2 ,C^, . . . ,C n parts of the program. To perform this task, it is necessary 
to keep all program parts organized for execution. This is done by using 
the proper push buttons on the Function Key Box (see Section 3 of Part II). 

To simplify the loading and executing of program parts, we use a 
command language that is executed interpretively ,• and which keeps the 
Cj's and their associated files organized for manipulation. To accomplish 
this kind of organization, a storage management scheme must be available 
that is general and easy to implement. The OLMS uses an algorithm called 
the T-Algorithm to accomplish this management scheme. This algorithm 
requires that the computer system being used to implement the interpreter 
have an overlay structure. The T-Algorithm is outlined below: 



T-Algorithm 

T1 : Create a simple overlay program that 

calls the interpreter for the conmand 
language. 

T2: The program created in step Tl, will 

be a part of every that we wish to 
interact with at execution time. Note, 
a call to the program created in step Tl 
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T3 : 



may appear anywhere in C. . 

To start the system, make a call 
to the program created in step Tl. 



The T-Algorithm functions so that any task executed from the interpreter 
will overlay the interpreter. This algorithm allows a user to develop 
a good simulation program under most existing operating systems. Finally, 
this algorithm generates flexibility for interaction with the Rand Tablet 
and Function Keys. 

The sections of this report that follow are an implementation of the T- 
Algorithm. To understand the details of how the T-Algorithm works, refer 
to the following discussion and diagrams. 



An Example of the T-Algorithm 

Let INTERP be the name of a routine that calls an interpreter (called 
INTER) . 

Let DMAIN be the name of a routine that calls INTERP to start the 
system. 

Let be a part of a program that displays data. Assume that INTER 
and the C . ' s are on disk. The diagrams below show how the memory would 
look under the T-Algorithm. 




DISC MEMORY 



The interpreter is now in memory and we can interact with it in any 
desired way. From the interpreter we will now allocate the C^’s. 




The program to display data is now in memory. When we wish to return 
to our interpreter, a call to INTERP is executed. 




INTER 



The interpreter is now in memory and our algorithm may continue. 
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SECTION II 



A Syntactic Description of A 
Modeling Command Language 



Language Description 

The language in which a language may be defined is termed a metalanguage and 
must be uniquely distinguishable from the language being described. To 
formalize the definition in the metalanguage, each definition is given the 
form of a statement or construct, which is analogous to a formula. However, 
to accomplish some unique features of such a specification, the operators 
define a mode of construction, or concatenation. The following symbols are 
employed in this metalanguage: 

X the object named X 



The syntactic part of this paper is divided into sections. Each section begins 
with a section number. Within each section there are subsections. All sub- 
sections are numbered with subsection numbers. An example of the structure of 
the syntactic part of this paper is given below. 



can be formed from 




or (the exclusive or) 

Z is to be repeated at least i times but not more 
than j times. When i is omitted, its value is to 
be assumed to be 1, and when j is absent, its 
value is assumed to be infinity. 



1 . Section 



1 . 1 Subsection 



1.1.1 Subsection 



2 . 



Section 



The last subsection of every section contains the general semantics of the 
section. 

1. Letter, Special Character, Null, Blank, Digits, Sign Numbers 

1.1. Letter 



) 



< Letter >: 



A | B | C | D 



|Z 



1.1.1. Semantics 

Letters do not have individual meaning. They are used for forming variables. 

1.2. Special Character 

<special character> : = + | - 1 , 

1.3. Null 

<null>: = A 

1.3.1. Semantics 

Lambda is a string of length zero. 

1.4. Blank 

<blank> : = 

1.4.1. Semantics 

Blank is a "special character" of length one. 

1 1.5. Digits 

<digit> : = 0|1|2|3|4|5|6|7|8|9 

1.5.1. Semantics 

Digits are used for forming numbers and variables. 

1.6. Sign 

<sign>: = <null> |+|- 

1.6.1. Semantics 

A null sign is an implicity plus sign. 

1.7. Exponent 

2 

<exponent>: = E<sign>{<digit >} i 

1.8. Numbers 

1.8.1. Integer 



10 

<integer>: = {<digit>} 1 
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1.8.2. Decimal Number 



<decimal number> : = {<digit> {<digit>} 1 0 “ n 

1.8.3. Fraction 

in 

<£raction>: = {.{<digit>} 1 } 

1.8.4. Number 

<Number> : = <fraction> | <decimal number> | 

<integer>|{< decimal number > |<fraction>>i 
l 

{< exponent;* }j 
1.1.0. Signed Number 

<Signed number> :=sign><number> 

2. Variable 

<Variable> : =<Name> 

3. Name 



4 

<Name> : = <letter> }<letter> | ^digit> } 0 



4. Save statement 

<Save statements = SAVE{<Blank > | , } i 

.1 l 

{ <Directory statement> } : { <Blank> | , } j 
l 

{<File>}i 

4.1. Semantics 

Files may be saved on magnetic tape, Dectape, and Disk. All files 
are saved from the system buffer. The system buffer is a named 
common block called COMBUF. The buffer size is three hundred and 
eighty four words. 

5. Get statement 

l 

<Get statement;* := GETF{<Blank>} 

1 . i i l 

{<Directory statement>}i {<Blank>|, }i(File>)i 

5.1. Semantics 

Files are read from magnetic tape, Dectape, and Disk into the system 
buffer . 
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- 6. Delete statement 

, l 

^Delete statement> : = DELE{<Blank > | , } \ 

l .1 l 

{<Directory statement^ {<Blank>|,} 1 {<File>} 1 

6.1. Semantics 

Files may be deleted from magnetic tape, Dectape and Disk. Any 
deleted can never be recovered. 

7. System Directory statement 

<System Directory statement : = NEWD 

. l 

{<Blank>|,)i <Device Name> 

7.1 Semantics 

The system directories can be written on magnetic tape, Dectape, and 
Disk. The file and library directory are called the system directories. 

8. Directory Input statement 

<Directory input statement : = DIRM 

, 1 1 i 1 11 

{<Blank>|,} 1 {<Device Name^ { { <Blank> | , } j { directory statement^! } 0 

8.1 Semantics 

If the device name is magnetic tape the magnetic tape directory will be 
read. A device name for Dectape and Disk will read the system directories. 
When the system directories are read the old contents of the directory 
tables are destroyed. For magnetic tape the <Directory statement> must be 
omitted. 

9. Clear statement 

<Clear statements = CLEA {<Blank>|,)i 
l 

{<Device name >}j 

9.1. Semantics 

The Dectape and magnetic tape directories can be cleared. The disk 
directory cannot be cleared. 

10. Rename statement 




I 1 1 

<Rcnamc statements : = RENA {<Blank>|,}i {directory statemcnt>} 1 

1 i , l 

{<Blank>j,} 1 {<File>} 1 {<Blank>|,}j 

, l l 

TO {<Blank>|,} 1 {<File>}j 

10.1. Semantics 



Files may be renamed on magnetic tape, Dectape and Disk. 

11. Description statement 

description statement : = DSCR 

i 1 1 

{<Blank>|,} 1 {directory statement^] 

l l 

{<Blank>|,} 1 {<File>}i 

11.1. Semantics 

File descriptions of fifty characters may be printed in the system 
directories. 

12. Plant statement 

<Plant statement : = PLAN {<Blank>|,} 1 

l 

{ directory statement } i 
{<Blank> | ,} 1 {<File>} 1 

12.1. Semantics 



File names may be planted in the system directories . 

13. Assign statement 

l 

<Assign statement : = ASSN {<Blank>|,} 1 

{device name> | {CHAN{<Blank> | , } j 

1 i 1 

<Device name>} 1 | {OLMS{<Blank>} ,}i 
l l 

<Device name>} 1 } 1 



13.1. Semantics 



All devices are assigned at execution time. These assignments are 
only methods to update system device tables. CHAN in the commands 
tells the system where user chain files are stored. OLMS in the 
commands tells the system where the On Line Modeling System is 
stored. 

14. Ueassign statement 

. l 

<Deassign statements* = DEAS{<Blank> | , }j 

11 

(<l)evice name> | (CHAN{<Blank> | , >i<Device name}* 



| (OLMS{<Blank> 



l i l 

, } x <Device name>}i}! 



14.1. Semantics 



Devices are deassigned at execution time. All deassignments are 
local, that is, they update system tables. CHAN in the command 
deassigns all CHAIN Files. OLMS in the command deassigns the 
system. 

15. Model statement 

<Model statement : = { <Blank> | , } } MODEL! <Blank> | , } j 

<Name>(<Blank> | , >i (A| S | G > j 

15.1. Semantics 



The A refers to an analysis model. S and G refers to a simulation 
and generation model respectively. 



16. Start statement 

< Start statement> : 



STARTC <Model 



l l 

statement>) i (<Core statement^ 



16.1. Semantics 

If the core statement is omitted the model will be run in 20k 
including the root segment and Block Data. The core size can 
be determined when the CHAIN file is generated. The loader 
returns the CHAIN size when a file is generated. If this size 
is greater than 20k the core statement should be used. 





17. 



Process statement 
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<Process statement>: = PROCE{<Blank> j , 

l l 

{<Model statement;*^ {<Core statement;* }q 

17.1. Semantics 

The Process statement will process the output from a simulation 
model. The core rules are the same as those in 16.1. Remember 
these syntactic rules are defined on tables which are updated 
at execution time; therefore, a simulation model or generation 
model can be run using this statement if the name has been 
given with the A property in the ENTER statement. 

18. Kill statement 

<Kill statement : = KILLM <Model statement 

18.1. Semantics 



KILL is local in that it removes the model name from the model 
name table. 

19. Enter statement 

l 

<Enter statements = ENTER{<Model statement 

19.1. Semantics 



v 

•i 



The ENTER statement is used to put model names in the model name 
tables . 

20. Display statement 

<Disp statement : = DISP 

20.1. Semantics 

This command will bring in the OLMGS. The OLMGS is operated inter- 
actively using the RAND Tablet. All the capabilities of the OLMGS 
is utilized to fit any users needs. 

When the OLMGS is loaded a user can use it interactively as described 
in the flow diagram in Figure 2. The flow diagram is designed to 
show the appearance of the CRT at various stages. 

The commands available in this program are listed below. 

<D0UT statement : = DOUT 

<Top statement : = T 

<Yes-No statement;*: = {Y|N}* 
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A carriage return is used to proceed in multipicture displays. 

To understand the capabilities of the OLMGS read reference [12] . When 
multipictures are displayed the system block transfer function is used 
to move blocks of data. A diagram in figure 3 shows how the data is 
moved from one buffer to another to accomplish the multipicture display 
task. Consider the example given in Figure 1. 



SYSBUF 



WORKSY 




Figure 3 



The first picture displayed is the one in SYSBUF. The dotted lines show what 
happens when the second picture is displayed. The lines show what happens 
when the third picture is displayed. 

In Figure 2 the user flags the boxes that appear on the CRT on the Rand Tablet 
to accomplished certain yes-no task. A flag is defined as a close of the 
stylus switch followed by a release. To input parameters when they are asked 
for by the system via the Rand Tablet a user should place the pen on the Rand 
Tablet depress it and move it to the right horizontally watching the CRT for 
the value box. When the value box is set to the desired value then raise 
the pen. 

21. Create statement 

<Create statements = CREA 

21.1. Semantics 

This command must be given before the Disp command in 20. When CREA is 
realized certain flags are set to prevent the user from doing anything 
but DISP. This is done to let the user think a fraction of a second 
before he gives the DISP command. The CREATE flags can be turned off 
without giving the DISP command. 
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22. DOUT statement 



<DOUT statement> : = DOUT 

22.1. Semantics 

The DOUT statement is used to return to the OLMS from the OLMGS. A 
DOUT given at the level of the OLMS will turn off the CREATE flags. 

23. Top statement 

<Top statement : = T 

23.1. Semantics 

This statement allows a return to the beginning of an interactive 
process. 

24. Core statement 

<Core statement : = 1 

24.1. Semantics 

The core statement tells the OLMS to allocate the maximum amount of 
core available. If no core is available the appropriate message 
will be typed and the system will return for another command. 

25. Transfer statement 



25.1. Semantics 

The transfer statement is defined on the system buffer and a work 
area. The work area is a named common block called WORKSY which 
is seven hundred words long. The OLMGS uses this area when multiple 
pictures are being displayed on one graph. If multiple pictures 
are to be displayed by the OLMGS, it is necessary to place the 
picture data in the buffers at the OLMS level. The OLMGS is de- 
fined on the first 256 words of system buffer. The first 128 words 
are the X- values and the second 128 words are the Y- values. The 
buffers in figure 1 shows how the transfer command works. 

Suppose we have done a GETF of a file of data. Now the file of 
data is in the system buffer. Let us call this file c^. Give 
a transfer to put 04 in the work buffer. 



<Transfer statements = TRANS {<Blank >|,} 1 
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SY°BUF ' JL28 words 


128 words 


fa fa /a /a 

t 1 / 1 / 2 / 2 




WORKSY 128 words 


128 words 128 words| 128 words 





Figure 1 



TRANS, 1 

Now do a GETF of a file called Give another transfer. 

TRANS, 2 

We have two pictures data sets stored in WORKSY. Give another GETF and 
call the file data a^. We have a,, a-, , as the first, second and third 
picture sets respectively. 

26. Generation statement 

1 

Generation statement> : = GENET {<Model statement>} 1 

26.1. Semantics 

Every model will require a specific kind of input. To meet this 
need the generation statement has been implemented to operate as 
the START and PROCESS statements. This command is discussed in 
detail in section three of this paper. 

27. File 

1 , 1 

<File>: = <Name>{{. <EXT>} 0 (<Blank>|,} 1 

1 1 

{<Device name>}i}i 

27.1. Semantics 

The extension may be omitted. 

28 . Device Name 

1 

<Device Name>: = MTA{<Unit num> } 1 \ DTA 

1 , 11 
{<Unit number>}i j {DSK{<Blank>} 1 }j 

28.1. Semantics 

These are the only devices on most systems that users may access. 

29. Property statement 

<Property statements = MrA|DTA|DSK 



O 
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9.1. Semantics 



The properties MTA, DTA, DSK are magnetic tape, Dectape and Disk 
respectively. 

30. Directory statement 

l 

<Directory statements = {L|F}j 

31.1. Semantics 

L and F refers to the system library directory respectively. 

31 . File Name 

l l 

<File Names = <Name>{ . <EXT>} 0 {<Blank> | ,}j 

] 

<Property statement>}i 

31.1. Semantics 

The extension may be omitted. 

32. List statement 

<List statements = LIST { <Blank> | , } i (A | S | M | D [ T | G | LA j FA} x j {L | F} i { <Blanlc> | , } i 
l 

<File Name>} 0 

32.1. Semantics 

This command will type out the system tables. 

The terminals are defined below. 



A List analysis models 

S List simulation models 

M List system mode 

D List devices assigned to system 

T List magnetic tape directory 

G List generation names 

LA List the whole system library directory 

FA List the whole system file directory 
L List a file name, extension, property and 
description from the system library 
directory 

F List a file name, extension, property and 
description from the system file directory 





33. Mode statement 
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<Mode statement> : = MODE{<Analysis statement> | <Simulation statement | 

<Function Key-light pen statements | <Expand statements | <Edit statements 

l 

<Null modes|<Build statements | <Generate statements } 1 

33.1. Semantics 

While designing the OLMS command language considerable thought 
was given to forcing the user to do what he wanted to do, but 
at the same time giving him considerable time to think and 
change his mind. To accomplish this task the language is de- 
signed with eight operation modes. In our implementation any 
mode can be entered from another mode. 

34. Analysis statement 

l 

<Analysis statements: = (<Blanks|,} 1 ANAL 

34.1. Semantics 

The commands available in this mode are listed below. 

<Mode statements: = MODE {< Analysis statements >. 

j 

| Simulation statements | <Function Key-light pen statements | <Expand statements 

iii i 1 

| <Edit statements | <Null modes |<Build statements | <Generate statements } 1 

<List statements: = LIST{<Blanks | , } x 

(A | S | M | D | T | G | LA | FA} \ | (<Blanks | ,}J 
l l 

(<File names} 

l 

<Save statements: = SAVE(<Blanks | , } 1 

l l 

{ <Directory statements } x { <Blanks | , } ± 

(<Files}J 

l 

<Get statements: = GETF{<Blanks | , } ^ 

l . l 

(<Directory statements } 1 {<Blank| , }x 
l 

(<Files } 1 

l 

<Delete statements; = DELE (<Blanks|,} 1 

1 i 1 > 

{ <Directory statements } 1 { <Blanks | , } x } 

{<Files } 1 

l 
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<Rename statement> : =RENA{<Blank> | , } j 

l l 

{ <Directory statement } i { <Blank> I , } 1 

1 i 1 i 1 1 i 1 

{<File>} 1 {<Blank> j , } 1 TO{<Blank> | , } j {<Directory statement>} 1 {<Blank> | , }j 
l 

{<File> >i 

<Description statement;- :=DSCR 

1 1 
{ <Blank> | , } 1 { <Directory statement> } j 

{ <Blank > | , } 1 { <File > } i 

<Help statements :=HELP 

<Disp statement> :=DISP 

<Create statement :=CREA 

<DOUT statement> :=DOUT 

l , l 

<Trans£er statement> :=TRANS{<Blank > | , } ± {112}! 

l 

<Process statement> :=PROCE{<Blank > | , } i 

l l 

{<Model statement>} 1 {<Core statement>} 



35. Simulation statement 

. l 

<Simulation statement :={<Blaak> | , } jSIMM 
35.1. Semantics 

The commands available in this mode are listed below. 

l 

<Start statement> := START {<Model statement}! 

l 

{<Core statement>} 0 

<Mode statement> :=MODE{<Analysis statement> | 

<Simulation statement> | <Function Key-light pen statement> | <Expand statement> 

l 

<Edit statement> | <Null mode>|<Build statement> | <Generate statement;- }} 

l i ,1 

<List statement> :=LIST{<Blank> | ,} 1 {A|S|M|D|T|G|LA|FA}i{<Blank> | ,>i 

l l 

{<File Name>} 1 } 0 

36. Build statement 
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<Build statement> :={<Blank> | , ijBUIL 
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36.1. Semantics 



When the build made is entered the system initilization flag is set to 
initialize the system in null mode when it is reloaded. The system 
exits to the computer monitor. 

37. Null Mode 

<Null mode>: = {<Blank>|,} 1 Null 

37.1. Semantics 

The commands available in this mode are listed below. 

<Mode statements = <MODE{<Analysis statement j Simulation statement |< 
Function Key-light pen statement | <Expand statement j <Edit statement> 

|<Null Mode>|<Build statement> | <Generate statement>} 1 

<List statement>: - LIST{<Blank> | , } x { A| S j M | D | T | G | LA | FA} 1 1 {L | F } 1 {<Blank> | , } 

l i 

{<File name>) ! }o 

, l l 

<System Directory statement> : = NEWD{<Blank> j , } 1 {<Device name>} 1 

<Directory input statement>: = DIRM 

. i l.l 

{<Blank>| ,} 1 {<Device name>} 1 {{<Blank> | 

l l 

{ <Directory statement> } i } q 

<Clear statement>: = CLEA{ <Blank> | , } j 

l 

{<Device name>} 1 

l 

<Plant statement> : = PLAN{<Blahk > | , } 1 

l 

{ <Directory statement > } i 

l l 

{ <Blank> | , } i { <File > } i 

<Assign statement>: = ASSN{<Blank> | , } i 

. l 

{ <Device name> | { CHAN{ <Blank > | , } ^ 

l . .1 

<Device name >} ^ \ {OLMS{<Blank > | , } i 

l l 

<Device name> } i } i 

<Help statement> : = HELP 

<End statement>: = ENDS 



38 . Expand Statement 
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<Expand statement : = {<Blank> | ,}iEXPA 

38.1. Semantics 

The commands available in this mode are listed below. 

<Help statements = HELP 

<Mode statements = MODE{< Analysis statement | <Simulation statement>|< 
Function Key-light pen statement | <Expand statement> | <Edit statement> 
|<Null mode>|<Build statement> | <Generate statement^ 

<List statements = LIST{<Blank> | , } x 
{A|S|M|D|T|G|LA|FA}J | (L J F } J { <Blank> | 

i l 

{<File name>} 1 } 0 

1 

<Enter statement> : = ENTER! <Medel statement>} 1 

i 

<Kill statements =KILLM {<Model statement 

39. Edit statement 

l 

<Edit statement> : = {<Blank> | ^EDIT 

39.1. Semantics 

Network models will be editted in this mode. Any variable 
in Block Data which is the memory of the OLMS can be changed. 

40. Generate statement 

. l 

<Generate statements = [<Blank> | ,>iGENN 

40.1. Semantics 

The commands available in this mode are listed below. 

<Help statements = HELP 

<LIST statements = LIST{<Blank> | ,}}{A|SIM|D|T|G|LA 

| FA} j | { L | F } i { <Blaiik> | ,}J{<File name>}J}o 

<Mode statement>: = MODE {< Analysis statement | <Simulation statement 

|<Function Key-light pen statement | <Expand statement | <Edit statement> 





a 





i i 1 

1 <Null mode> j <Build statements J <Generate statement^ j 

l 

Generation statement> : GENET {<Model statement>} 1 


41. 


Function Key-light Pen statement 

l 

<Function Key-light pen statement> :={<Blank> | , } iFKLP 


41.1. 


Semantics 

In this mode certain commands in the command language are 
assigned to function keys. This mode is discussed in Section III 


42. 


Unit Number 

<Unit Number>:= 0|l|2|3|4|5|6|7 


43. 


Unit Num 

<Unit Num> : = 0 1 1 j 2 


44. 


EXT (File extension) 

, . 2 
<EXT> :=<Null> | <Letter>{<Letter> | <digit>) o 


45. 


Help statement 
<Help statement :=HELP 


45.1. 


Semantics 

Help is a program that teaches a user how to use the OLMS. 
The commands available in this program are listed below. 

<Top statement :=T 

<Yes-No statement> :={Y|N}i 

<Dout statement> :=DOUT 


46. 


END statement 

<END statement :=ENDS 


46.1. 

47. 


Semantics 

This command terminates the system. 

Yes -No statement 

<Yes-No statement> :={Y|N} ! 
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SECTION III 



GENERATION OF INPUT 

The Rand Tablet and Function Keys are excellent input devices for interactive 
work with interactive systems. To use these devices economically it is necessary 
to have a flexible system. This flexibility is obtained in the OLMS by allowing 
any type of INPUT program using the Rand Tablet and Function Keys. 

The OLMS has a standard graphical input using the Rand Tablet, The flow chart 
in figure 1, shows the sequence for inputting a graph and collecting the graphical 
data input. The large blocks show how the CRT looks at certain levels in the 
program. 

Commands available are listed below. 

<Top command> : =T 
<Dout command> : =DOUT 



A carriage return will return to the data bounds input level. 

This implementation allows a user to collect 128 X-values and 128 Y-values. 

FUNCTION KEY LIGHT PEN MODE 

In this mode various commands of the Command Language are assigned to 

Function Keys. The assignments are listed below. 

COmAND KEY ASSIGNMENT 

SAVE 1 

MODE 2 

DISP • 3 

GETF 4 

DELE 5 

TRANS 6 

START 7 

CLEA 8 

RENA 9 



<Transfer command> :=TRANS{<Blank> | ,}^{1|2}^ 




COMMAND 



KEY ASSIGNMENT 



DSCR 10 
PLAN 11 
ASSN 12 
DEAS 13 
PROCE ' 14 
GENET 15 
DOUT 16 
KILLM 16 
ENTER 16 



This mode is operated by depressing a key associated with a particular 
command and following the program sequence. Whenever a decision is needed 
a small box will appear adjacent to the statement on the CRT. The Rand Tablet 
is used to flag these boxes. If data input is requested on the teletype the 
CONTINUE will be typed before the system is ready for input. The program 
sequence is listed in figure 2. The large blocks show the appearance of the 
CRT at a particular level in the program. 



TOP 




SET MARGIN 



Iyes 1 

Ino 1 



LEFT MARGIN 
RIGHT MARGIN 
BOTTOM MARGIN 
TOP MARGIN 



VALUE 

BOX 



SET GRAPH DENSITY 



VALUE 

BOX 





GRAPH TYPE 
LINEAR-LINEAR 

LOG -LINEAR 

LINEAR-LOG 

LOG -LOG 







.> 



•' > 



) 
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FUNCTION 


KEY MODE 


COMMAND 


KEY 


SAVE 


1 


MODE 


2 


DISP 


3 


GETF 


4 


DELE 


5 


TRANS 


6 


START 


7 


CLEA 


8 


RENA 


9 


DSCR 


10 


PLAN 


11 


ASSN 


12 


DEAS 


13 


PROCE 


14 


GENET 


15 


DOUT 


16 


KILLM 


16 


ENTER 


16 


. 


£ 



<2 — ( depress m) 
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SECTION IV 




4.1 Example 

This example shows how the OLMS can be used to study the input 
output relations of the Hodgkin-Huxley Nerve Equations. The 
differential equations are solved using the Predictor -Corrector 
Method with the Runge-Kutta Method as a starter. 

The equations are given below under the heading Hodgkin-Huxley 
Equations. Any details about these equations can be found in 
references [29] and [28]. 

Hodgkin-Huxley Equations 

(1) I=CV+g Na m 3 h (V+115) +g K n 4 (V- 12) +g L (V+10 . 5989) 

(2) m =<J> [ (l-m)a m (V) -m$ m (V) ] , 

(3) h =<(>[(l-h)a h (V)-h3 h CV)] , 

(4) n =4[(l-n)a n (V)-nS n (V)]. 
where 

a m (V) =0 . 1 CV+2S3 [exp(^i) - 1 ] ' 1 , 

6 m (V)=4exp(V/18) , 

a h (V)=0.07exp(V/20) , 

e h (V) = [exp(^) + l]- 1 , 

a n (V) =0 . 0 1 (V+10 j [exp - 1 ] ‘ 1 , 

B n (V)=0.125exp(V/80) , and 
^ =3 (T-6.3)/10 

The Predictor-Corrector used is given as follows; 
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Predictor 



Given m points 

x ,x , ,x 

n-m+1 ’ n-1’ n 

£ ,£ , ,£ and Newton's Integrating Newton's backward 

n-m+1 n-1 n 

interpolation formula we have: 

m-1 

y , , = y +h Z B . A-^f 
n+ l n j=o J 

where 

x-x 

z variable at x is z= —r — and 
n h 

j-1,2, ,m-l. 



Corrector 

Given: 

x ,x ,x , , 

n-m+1’ ’ n’ n+l 

f m+1’ ,£ ,f ^ and Newton's Integrating Newton's backward 

interpolation formula we have: 

m-1 * . 

VrV h . z n V Wj 

J=0 J J 

where 

A "l 0 ^ 

B 0 =1 , B C“ J /_i C j)dz, 

x-x 

z variable at x is z= -t- — and 
n h 

j=l,2, ,m-l 



This model is given in algorithmic form and a program listing is given 
for the solution of the model. This model can be extended by modifying 
the algorithm and using two dimensional arrays to solve the partical 
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differential equations in reference [28]. Someone may claim the 
algorithm uses too much storage. An algorithm can be developed 
that checks for vector overflow. If the vectors are full the data 
in the vectors from the beginning up to the total length minus three 
can be outputted. The three data points can be moved to the 
beginning of the vectors. By adjusting the appropriate pointers the 
algorithm will think it is starting. It is not necessary to move 
the data if linked data structures are used. The stepsize can be 
changed by creating an algorithm that starts the Runge-Kutta when 
some condition occurs, perhaps when the Predictor-Corrector fails to 
converge. This can be done by backing up one step and starting the 
Runge-Kutta with the new stepsize for some number of points greater 
than or equal to three in our case. Adjust the appropriate pointers 
and give control to the Predictor-Corrector. 

The Model in Algorithmic For m 

In thealgorithms that follow it is assumed that the relationship of 
functions and subroutines exist. There are two classes of operators 
called Operator and Parameter and Transfer Operator. The Operator and 
Parameter class is used to give input to the algorithms . The Transfer 
Operator is used to transfer data to the OLMS system buffer. All 
operators are defined below: 

Operator and Parameter Class. 

CAP - membrane capacitance 

UCON - turn off all variables held constant 
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DX 



« 



DX - 
COLL - 


distance steps ize W 

data collection steps ize used to transfer data to the OLMS 
system buffer. 


RUN - 


• start algorithm 


TOWM - 


• T 
m 


TOWH - 


T h 


TOWN - 


• T 

n 


TEMP - 


■ temperature (°C.) 


M 


■ sodium activation 


H 


• sodium inactivation 


N 


• potassium activation 


T 


• initial time 


V 


- initial voltage 


CUR ■ 


■ stimulating current pulse ) 


DUR ■ 


- duration of current pulse 


SZ 


- integration steps ize 


EPS • 


- convergence term 


STAR ■ 


- number of starter values 


FT ■ 


- number of solution points desired 


MITE • 


- maximum number of corrector iterations per step 


MO 


- hold sodium activation constant 


HC 


- hold sodium inactivation constant 


NC 


- hold potassium activation constant 
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Transfer Operator Class 
TV - time .vs. voltage 
TM - time .vs. sodium activation 
TH - time .vs. sodium inactivation 

TN - time .vs. potassium activation 

VM - voltage .vs. sodium activation 

VH - voltage .vs. sodium inactivation 

VN - voltage .vs. potassium activation 

DI - distance .vs. current 

TI - time .vs. current 

The global variables in the algorithms are defined as follows 
H - integration stepsize 
EPS - convergence term 
TEMP - temperature 

NC - total number of solution points desired 
N - number of starter values 

MAX - maximum number of corrector iterations per step 

PHI - <)> for sodium activation 

PHH - <j> for sodium inactivation 

PHN - <j> for potassium activation 

ICO - counter 

DX distance stepsize 

ARAD - radius of axon 

RR - specific resistance of axoplasm 

SIS - amplitude of the stimulating current 

CURR - current 



DIST 



- distance along axon 
DURRAT - duration of current pulse 
CAP - membrane capacitance 
CHAR - control vector 
YHH - vector os sodium inactivation 

YMM - vector of sodium activation 

YNN - vector of potassium activation 

X - vector of time 

SYSBUF - OLMS system buffer 

The R-Algorithm is used to input data and transfer data to the 
OLMS system buffer. The A-Algorithm is the main algorithm that uses 
the V,FUM,FUH > FUN,P,C, and FN algorithms to integrate the equations. 



R-Algorithm 
Subroutine RUNMOD 



Rl. [Give Operator and Parameter] 
VAR •'-INPUT 
VALUE ^INPUT 



R2. [What is the Operator?] 

if VAR.EQ. 'CAP' then CAP •'-VALUE , go to Rl, 
Otherwise 




if VAR.EQ. 'UCON' then CHAR(22) 0,CHAR(23) <- 0,CHAR(24) 0, 

go to Rl, 

Otherwise 

if VAR.EQ. 'DX' then DX <- VALUE, go to Rl, 

Otherwise 

if VAR.EQ. 'COLL' then VALCOL VALUE, go to Rl. 

Otherwise 

KCT +■ 0 (initialize control count) 

LOOP: KG 1 to 20 (loop to set KCT) 

KCT •<- KCT+1 

if VAR.EQ. CHAR (KCT) then go to LO, 

Otherwise go to LOOP 
Operator Error, go to Rl. 

LO: if KCT.EQ.l then go to R3, 

Otherwise 

if KCT. EQ. 2 then TOWM VALUE 
if KCT.EQ.3 then TOWH «- VALUE 

4-7 
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A-Algorithm 
Subroutine ADAMS 



« 



Al. [Generate Starting Points] 

N +■ The number o£ starting values 

MAX *■ The maximum number o£ iterations per step 

EPS ■*- The convergence region 

H «- Stepsize 

NC *■ Maximum number o£ steps 

NCMN «- NC-N The number of steps to use the predictor-corrector 
I «- 0 Initialize the predictor-corrector step counter 
CALL RUNFIT Generate N starting points 
ICO <- 0 

LOOP: ICO <- ICO+1 

if ICO <_ N then go to LOP 
Otherwise ICO -<- ICO-1, go to LO. 

LOP: F(ICO) VOLT (Y( ICO)) 

FMM(ICO) FUM(YMM(ICO) ) 

FHH(ICO) «- FUH(YHH(ICO)) 

FNN(ICO) «- FUN(YNN(ICO) ) 

X(IC0+1)=X(IC0)+H 
LO: if ICO < 3 then ERROR. 

A2. [Start the predictor- corrector] 

ICO «- N+I 

Initialize convergence switches 
IVSW 0 , IMSW 0 ,IHSW ^ 0 , INSW ^ 0 
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R7. [Transfer Time .VS. Sodium Activation] 

Same as R6 after replacing Y(I) in A2 by YMM(I) and L00P1 by L00P2. 

R8. [Transfer Time .VS. Sodium Inactivation] 

Same as R6 after replacing Y(I) in A2 by YHH(I) and L00P1 by LOOPS. 

R9. [Transfer Time .VS .Potassium Activation] 

Same as R6 after replacing Y(I) in A2 by YNN(I) and L00P1 by L00P4. 

RIO. [Transfer Voltage .VS. Sodium Activation] 

Same as R6 after replacing X(I) in A1 by Y(I), Y(I) in A2 by 
YMM(I) and L00P1 by LOOP5. 

Rll. [Transfer Voltage .VS. Sodium Inactivation] 

Same as R6 after replacing X(I) in A1 by Y(I), Y(I) in A2 by YHH(I) 
and LOOP1 by LOOP 6. 

R12. [Transfer Voltage .VS. Potassium Activation] 

Same as R6 after replacing X(I) in A1 by Y(I), Y(I) in A2 by YNN(I) 
and LOOP1 by LOOP 7. 

R13. [Transfer Distance .VS. Current] 

Same as R6 after replacing X(I) in A1 by DIST(I), Y(I) in A2 by CURR(I) 
and LOOP1 by L00P8. 

R14. [Transfer Time. VS. Cur rent] 

Same as R6 after replacing Y(I) by CURR(I) and L00P1 by L00P9. 

R15. [Return to the OLMS] 



-r, 



o 




3c? 



R4. [Transfer data to the OLMS system buffer] 



VAR <- INPUT 

R5. [Check Transfer Operator] 
VALTEM «- VALCOL 

J 1 



if 


VAR.EQ. 


'TV' 


then 


go 


to 


R6 


if 


VAR.EQ. 


'TM' 


then 


go 


to 


R7 


if 


VAR.EQ. 


'TH' 


then 


go 


to 


R8 


if 


VAR.EQ. 


'TO' 


then 


go 


to 


R9 


if 


VAR.EQ. 


'VM' 


then 


go 


to 


RIO 


if 


VAR.EQ. 


'VH' 


then 


go 


to 


Rll 


if 


VAR.EQ. 


’VN* 


then 


go 


to 


R12 


if 


VAR.EQ. 


*DI * 


then 


go 


to 


R13 


if 


VAR.EQ. 


•TI * 


then 


go 


to 


R14 



Otherwise 



Transfer Operator Error, go to R4 
R6. [Transfer Time. VS. Voltage] 

L00P1: Increment I by 1, if 1=500 then finished otherwise, 
if J.GT.128 then go to LOOP1 
if X(I).NE. VALTEM then go to LOOP1 

Otherwise 
SYSBUF (J) +- X(I) 

SYSBUF(J+128) «- Y(I) 

J -h J+l 

VALTEM -h VALTEM+VALCOL 

go to L00P1, go to R15 (after the loop is complete) 
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A2. 
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if KCT.EQ.4 
if KCT.EQ.5 
if KCT.EQ.6 
if KCT.EQ.7 
if KCT.EQ.8 
if KCT.EQ.9 
if KCT.EQ.10 
if KCT.EQ.ll 
if KCT.EQ.12 
if KCT.EQ.13 
if KCT.EQ.14 
if KCT.EQ.15 
if KCT.EQ.16 
if KCT.EQ.17 
if KCT.EQ.18 
if KCT.EQ.19 
if KCT.EQ. 20 
go to Rl. 

R3. [Start Model] 



then TOWN <- VALUE 
then YMM(l) - VALUE 
then YHH(l) <- VALUE 
then YHH(l) VALUE 
then YNN(l) VALUE 
then X(l) +■ VALUE 
then Y(l) «- VALUE 
then SIS «- VALUE 
then DURRAT «- VALUE 
then H VALUE 
then EPS ■«- VALUE 
then N VALUE 
then NC VALUE 
then MAX VALUE 
then CHAR (2 2) «- CHAR(21) 
then CHAR(23) CHAR(21) 
then CHAR(24) •*- CHAR(21) 



PHI DEXP (.1098612289D0* (TEMP-6. 3DO)) 
PHH + PHI 
PHN PHI 

if TOWM.NE. 1 then PHI + PHI/TOWM 
if TOWH.NE. 1 then PHH + PHH/TOWH 



if T0WN.NE.1 then PHN PHN/TOWN 
CALL ADAMS (start integration) 



Check for constant conditions 



if CHAR(22) .NE.0.0 then IMSW <- 1. 
if CHAR(23) .NE.0.0 then IHMS <- 1 
if CHAR (24) .NE.0.0 then INSW <- 1 
Initialize corrector iteration counter 
M -<- 0 
NPI <- N+I 

A3. [Compute f N+J = f (x N+I ,Y N+I ) ,etc. ] 

F (N+I) «- VOLT (Y (N+I)) 

Cehck for constant conditions 
if CHAR(22) .NE.0.0 then go to B 
Otherwise, ICO «- ICO-1 

FMM(N+I) «- FUM(YMM (N+I)) 

ICO + ICO+1 

Check for constant conditions 
0: if CHAR(23) ,NE. 0 . 0 then go to G 
Otherwise, ICO *■ ICO-1 

FHH(N+I) «- FUH(YHH (N+I) ) 

ICO «- ICO+1 

Check for constant conditions 
L: if CHAR(24) .NE.0.0 then go to LI 
ICO «- ICO-1 

FNN(N+I) «- FUN(YNN (N+I)) 



ICO i- ICO+1 



h 






A4. [Compute predicted value for Step I] 

if IVSW.NE. 1 then Y(NPI+1) <- PREDIC(Y,NPI,H,F) 
if IMSW.NE. 1 then YMM(NPI+1) <- PREDIC (YMM , NPI , H , FMM) 
if IHSW.NE. 1 then YHH(NPI+1) «- PREDIC (YHH, NPI ,H,FHH) 
if INSW.NE. 1 then YNN(NPI+1) «- PREDIC (YNN, NPI, H,FUNN) 
T: if IVSW.EW.O then go to E 
Otherwise , 

A5. [Set up for corrector] 



if IMSW.EW. 0 then FMM(NPI+1) 
if IHSW.EQ.O then FHH (NPI+1) 
if INSW.EQ. 0 then FNN (NPI+1) 



FUM(YMM (NPI+1) ) 
FUH (YHH (NPI+1)) 
+• FUN (YNN (NPI+1)) 



A6. [Compute corrected solution value for step I, iterati on M] 
if IVSW.EQ. 0 then YVC +■ CORREC(Y,NPI ,H,F) 
if IMSW.EQ. 0 then YC +• CORREC (YMM, NPI ,H,FMM) 
if IHSW.EQ.O then YCH <- CORREC ( YHH, NP I ,H,FHH) 
if INSW.EQ. 0 then YCN «- CORREC (YNN, NPI ,H,FNN) 

A7. [Test for convergence] 

DELV +- |YVC-Y(NPI+1) | 

DELY |YC-YMM(NPI+1) | 

DELH |YCH-YHH(NPI+1) | 

DELN | YCN- YNN (NPI+1) | 

if IVSW.EQ. 1 then go to 01 
if DELV-EPS <_ 0 then go to 02 
if DELV-EPS > 0 then go to 01 






02: Y(NPI+1)“- YVC ; good value converges 
IVSW 1 

01: i£ IMSW.EQ.l then go to 03 
if DELY-EPS<_ 0 then go to B1 
if DELY-EPS > 0 go to B2 
Bl: YMM(NPI+1) «- YC ; good value converges 
IMSW «- 1 

B2: if IHSW.EQ.l then go to G1 
if DELH-EPS £ 0 then go to G2 

if DELH-EPS > 0 then go to G1 

G2: YHH(NPI+1) «- YCH ; good value converges 
IHSW 1 

Gl: if INSW.EQ. 1 then go' to T1 
if DELN-EPS £ 0 then go to T2 

if DELN-EPS > 0 then go to T1 

T2: YNN(NPI+1) «- YCN ; good value converges 
I NSW «- 1 

Tl: if M-MAX £ 0 then go to El 
if M-MAX > 0 then go to E2 
El: if IVSW.NE. 1 then Y(NPI+1) «- YVC 
if IMSW.NE. 1 then YMM(NPI+1) <- YC 
if IHSW.NE. 1 then YHH(NPI+1) <- YCH 
if INSW.NE. 1 then YNN(NPI+1) <- YCN 
if IMSW+IHSW+INSIV+IVSW.EQ. 4 then go to S 
Otherwise, 

M *♦* M+l 
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go to T 



E: ICO +■ ICO+1 



IMSW <- 0 
1HSW +- 0 
INSW +- 0 

F(NPI+1) +- VOLT(V(NPI+l)) 

ICO +• ICO- 1 
go to A5 

B: [Constant YMM] 

YMM( ICO+1) +• YMM (ICO) 
FMM(IC0+1) FMM(ICO) 
go to 0 

S: if I-NCMN ^0 then go to SI 

if I-NCMN > 0 then go to S2 

SI: I +• 1+1 

X(NPI+1)=X(NPI)+H 
go to A2 

G: [Constant YHH] 

FHH( ICO+1) +- FHH(ICO) 

YHH (ICO+1) +- YHH (ICO) 
go to L 

LI: [Constant YNN] 

FMM( ICO+1) +- FNN(ICO) 
YNN(IC0+1) +- YNN (ICO) 
go to S 

S2: Terminate algorithm. 

E2: NC +- 1 




Fails to converge. 



FN- Algorithm 
Function FN (X) 



FN1. [Initialize BCOEFF vector] 
BCOEFF(l) + -.083333333 
BCOEFF (2) 4 - .0166666667 

BCOEFF (3) .0238095238 

BCOEFF (4) + .025 

BCOEFF (5) .025252525 

BCOEFF ( 6 ) <- .02531136 

BCOEFF (7) *► .02533 

FN2. if |X|-1. < 0 then go to FN5 
if jX | - 1 . ^0 then go to FN3 
FN3.FN ■«- X/ (EXP (X) - 1 . ) 

FN4. Terminate 
FN5.X2 «- -X*X 

IOFLO ♦ OVERFL (IOFLO) 

FN 4- 1. 

A 4- - ,5*X 

IOFLO 4- OVERFL (IOFLO) 
if IOFLO.EQ.l then terminate 
FN6.FN l.+A 
FN7 . [LOOP] 

LOOP: 1=1 to 7 

A 4- BCOEFF (I) *X2*A 
IOFLO 4- OVERFL (IOFLO) 



f 



i£ IOFLO.EQ.l then terminate 
FN «• FN+A 
go to LOOP 

terminate (after loop is complete) 

This is not the only way to implement this model using the OLMS 
however, it is one of the simplest. 




0 



V- Algorithm 



Function Volt (Y) 

V1.T0LT «- -120 .*YMM(IC0)**3*YHH(IC0)* (Y-115 .) 

-36 . *YNN(ICO) **4* (Y+12 . ) - .3* (Y-10. 598921) +CURRQCO) 
V2. if CAP.NE.l . then TOLT <- TOLT/CAP 
V3. VOLT <- TOLT 
V4. Terminate 
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FUH- Algor i thm 



Function FUH(YHH) 

FUH1.I00 ICO+1 

FUI-I2 . FUH 4 . PHH* ( (1 . -YHH) * . 07*EXP (- . 0 5*Y ( ICO) ) -YHH/ ( 1 . •'■EXP (3 . - . 1*Y ( ICO) ) ) ) 
FUH3.ICO •». ICO-1 



FUH4. Terminate 



FUM-Algorithm 



Function FUM(YMM) 

FUM1.IC0 •*- ICO+1 

FUM2.FUM <- PHI* ( (1 . -YMM) *FN (2.5-. 1*Y (ICO) ) -YM4*4 . *EXP (-Y (ICO) /18 . ) ) 
FUM3.ICO <- ICO-1 
FUM4 . Terminate 
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FUN- Algor ithm 



* 



Function FUN(YNN) 



FUN 1. 1 CO «- ICO+1 

FUN2.FUN -<- PHN*((1.-YNN)*.1*FN(1.-.1Y(IOO))-YNN*.125*EXP(-.0125*Y(ICO )) j 
FUN3.IC0 <- ICO-1 
FUN4. Terminate 



r 

r 

r 
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P -Algorithm 



Function PREDIC(Y,NP1 ,H,F) 

Pl.PREDIC Y (NPI) + (H/ 24 . ) * ( 55 . *F (NPI ) - 59 . *F (NPI - 1) +37 . * 

F(NPI-2) -9 .*F (NPI-3) ) 

P2. Terminate 






C- Algorithm 



Function CORREC(Y,NPI ,M,F) 

C1.C0RREC «- Y (NPI) + (H/ 24 . ) * (9 • *F (NPI+ 1) +19 . *F (NPI ) - 5 . *F 
(NPI - 1) +F (NPI - 2) ) 

C2. Terminate 



I 
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FN- Algor ithm 



Function FN(X) 

FN1 . 

FN2. if | X | - 1 . < 0 then go to FN3 
if | X | - . 1 0 then go to FN5 

FN3.FN f X/ (EXP (X) - 1 . ) 

FN4. Terminate 
FN6.X2 -X*X 

IOFLO OVERFL(IOFLO) 

FN7.FN «- 1. 

FN8.A «- - . 5*X 

This is not the only way to implement this model using the 
OLMS; however, it is one of the simplest. Several graphs are 
given to show the input output relationships of this model . 



RESULTS OP SIMULATION 

These graphs are given to show how the results of a model can be displayed 
using the OLMS. 




Action Potential 




Sodium Activation 






Potassium Activation 




Sodium Inactivation 
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Sodium Activation v.s. Voltage 




Sodium Inactivation v.s. Voltage 



O 

ERIC 



4-73 




Potassium Activation v.s. Voltage 
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The Division of Computer Research and Technology, NIH, will 
issue on an irregular basis technical documents which we believe 
will be of particular interest to the biomedical community. 

These reports will include detailed descriptions of relevant 
computer programs and instructions in their use (as well as some 
theoretical background) , in hopes that intereste ! scientists will 
be encouraged to gain firs.t-hand experience in applying them. 

In some cases, such reports may serve as foci around which DCRT 
will structure training courses to expand the knowledge and 
experience of NIH staff in applying computer science to problems 
of research and management. Circulation of these reports within 
the biomedical community broadly is, of course, encouraged. 






^ It). W.z 5. 

A. W. Pratt, M.D., Director, DCRT 



3 

ERIC 



79 



I 



U.S. DEPARTMENT OF HEALTH. 
EDUCATION & WELFARE 
OFFICEOF EDUCATION 
THIS DOCUMENT HAS REEN REPRO- 
DUCED EXACTLY AS RECEIVED FROM 
THE PERSON OR ORGANIZATION ORIG- 
INATING IT. POINTS OF VIEW OR OPIN- 
IONS STATED DO NOT NECESSARILY 
REPRESENT OFFICIAL OFFICE OF EDU- 
CATION POSITION OR POLICY 



THE ON LINE MODELING SYSTEM 



Part II, Programmers f Reference Manual 

by 

Edward Hill, Jr. 

Laboratory of Applied Studies 






o 




Division of Computer Research and Technology 
National Institutes of Health, Public Health Service 
Department of Health Education and Welfare, Bethesda, Maryland 20014 



I 



PROGRAMMERS ' 



REFERENCE MANUAL 



Section 



Table of Contents 



Page 



I INTRODUCTION 1-1 

II BASIC OLMGS PROGRAMMING ON THE DEC 340 DISPLAY 2-1 

SCALE 2-1 

INTENSITY 2-1 

MODE 2-1 

DATA WORD FORMATS 2-2 

PARAMETER MODE 2-2 

POINT MODE 2-2 

RASTER MODE 2-2 

CHARACTER MODE 2-3 

VECTOR MODE 2-3 

VECTOR CONTINUE MODE 2-3 

INCREMENT MODE 2-4 

SUBROUTINE MODE 2-4 

SYMBOL CONVENTIONS 2-4 

III OLMGS CONTROL ROUTINES 3-1 

SETTING ERROR POINTERS 3-1 

CHAIN FILE COMPATABILITY 3-1 

INTENSITY SELECTION: BRITEV, FATNTV, ITSCV 3-1 

SCALE SIZE: BIGV, SMALLV, CHARV 3-2 

IV GRAPHS 4-1 

Graphing Data: KWKPLT, PDP340 4-1 

V GRIDS 5-1 



Generating a Grid: GRID1V 5-1 

Examples of GRID1V Usage 5-3 

Computation of GRID1V Arguments: DXDYV 5-13 




i 



R ! 



*- 



Section 



Page 



V 

(cont) GRID1V Controls 5-14 

Grid Margin Variation: SETMIV, SETMOV 5-15 

Examples 5-16 

Providing for Special Label Characters: SETCIV, 

SETCOV 5-16 

Holding Margins from Graph to Graph: HOLD IV, 

HOLDOV 5-18 

Operational Details of GRID1V 5-19 

Determining Grid Boundary Positions 5-25 

GRID1V Error Procedure 5-26 

Examples of GRID1V Characteristics 5-27 

Example 1 5-27 

Example 2 5-27 

Example 3 5-32 

Example 4 5-32 

Example 5 5-32 

Log and Semi -Log Grids 5-35 

Restrictions on Logarithmic Mode 5-36 

Examples 5-36 

Building Special Grids: LINRV, NONLNV 5-39 

Axis Lines: XAXISV, YAXISV 5-44 

VI SCALING AND CONVERSION 6-1 

Basic Scaling Subprograms: XSCALV, YSCALV 6-1 > 

Example 6-2 

Basic Scaling Equations 6-2 

Conversion of Data: NXV, NYV, IXV, IYV 6-3 

Basic Conversion Equations 6-4 

Inverse Conversion 6-6 

Retrieval of Scaling Information: SCLSAV 6-7 

Resetting Scaling Information: RESCLV 6-8 

Nonlinear Scaling and Conversion 6-8 

Linear -Nonlinear Scale Mode Indicators: SMXYV, 

MSXYV 6-8 

Nonlinear Transformation: XMODV, YMODV 6-9 

Off-Scale Error Detection 6-10 

Set Conversion Error Indicators: SCERRV 6-10 

Example 1 6-11 

Example 2 6-11 

Saving and Resetting Error Indicator Cells: 

SERSAV, SERREV 6-12 







ii 



82 



Section Page 

VII PLOTTING 7-1 

Plotting Data 7-1 

Point Plotting Subprograms 7-1 

Plotting an Array: APLOTV 7-1 

Plotting Individual Data Points: POINTV 7-7 

Basic Plot-Print Subprogram: PLOTV 7-8 

Line Generation: LINEV 7-11 

VIII TITLING AND LABELING 8-1 

Titling and Labeling Subprograms 8-1 

Horizontal Printing: PRINTV 8-1 

Vertical Titles: APRNTV 8-2 

Fixed Point Labels: LABLV 8-4 

Number Conversion: NXV,NYV 8-7 

IX PRINTING 9-1 

OLMGS 340 PRINTING: TX340 9-1 

X MISCELLANEOUS ROUTINES 10-1 

Unpacking Characters: UNPACK 10-1 

Setting System Flags: FLAGS 10-1 

\ Drawing Vectors : PLTW5 10-1 

Display Premitives 10-2 



3 

ERIC 



iii 



83 



SECTION I 



INTRODUCTION 



To model anything it is necessary to analyze the results given by the 
model. In order to model it is important to have flexibility in both 
the modeling system and the modeling graphical subsystem. During the 
development of an On-Line Modeling System (Part I) , I found it desirable 
to develop an On-Line Modeling Graphical System (OLMGS) . This OLMGS 
is completely flexible with respect to graphical display. This flexi- 
bility will become obvious as this report is read. 

The most obvious application of the OLMGS is the rapid production of 
labeled graphs. Results in graphical form are usually much easier to 
analyze then are the same results in printed tabular form. An option 
that is available is the production of a picture that combines tabular 
data with a plot of the related curve. Still another use of the OLMGS 
techniques is the creation of diagrams and line drawings. A camera may 
be operated with this system under program control. 

At first glance, the list of the OLMGS routines on the following pages 
may seem formidable indeed, but the user should not feel baffled by 
the great volume of details. Some routines were developed for special 
purposes and have limited usage. Many are used principally as lower 
modules for more general subroutines, and will seldom be called 
directly by the programmer. 

This Programmer's Reference Manual has been organized in major sections. 
These are arranged according to the function of the routines described 
therein. The major sections are: Introduction, Basic OLMGS Programming 

on the DEC 340 Display, Control Routines, Grid Drawing Routines, 

Scaling Routines, Point-Plotting and Line Drawing Routines, Titling 
and Labeling Routines, and a Higher Level Printing routine. Information 
about the interpretation of the system subroutines has been included 
in the control section. 

The OLMGS is a set of Fortran IV and Macro- 10 Assembly Language routines 
that drive display premitives which were written by Mr. Harry Lewis. 

The Light Pen and Rand Tablet can be used in the manner described in 
Mr. Lewis' Reports, (see reference [3] and [5]). It is possible to 
implement this system on any computer if the rasters on the CRT are the 
same as the ones used in this report. If the rasters are different a 
few constants and uppt limits for loops can be changed to make the 
system compatible. 




ERIC 



1-1 

84 



SECTION II 



BASIC OLMGS PROGRAMMING 
ON THE DEC 340 DISPLAY 



The DEC 340 Display plots images on a grid 1024- horizontal by 1024- vertical 
points. There are therefore, over 1 million addressable positions. The 
center of any character may be plotted at any of these positions. The origin 
(0,0) is in the lower left comer, with X and Y increasing to the right and 
upward, respectively as shown in the figure below. The electron beam does 
not scan the face of the tube; its position is determined by the contents of 
the X and Y registers of the display. (If a data word specifies a point 
outside of the raster an edge violation occurs and stops the display) . 



SCALE 



0,1023 




1023,1023 



1023,0 



The scale setting determines the number of positions each succeeding spot is 
moved before it is intensified. It effects both the size and appearance of 
lines or symbols drawn in the vector, vector continue, increment, or character 
modes. At scale setting 11-, each point can be clearly distinguished. At 
scale setting 00„, lines ana symbols appear to be continuous. The point 
spacing is illustrated in the following table. 



INTENSITY 



Scale 


Point Spacing 


I tens i tv 


00 




Every 


01 


• o*o#o»o# 


2nd 


10 


• 

o 

o 

o 

• 

o 

o 

o 
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4th 


11 


• 

o 

o 

o 

o 

o 

o 

o 

• 


8th 



There are eight intensity levels available on the display, ranging from 00 ? 
which is barely visible, to 111., which is very bright. Note that scale ani 
intensity settings are interrelated. For example, if characters are drawn 
(with the character generator) at the lowest scale setting, and too high an 
intensity is used, they will be badly blurred. On the other hand, if many 
characters are to be displayed simultaneously or if the light pen is to be 
used, it is best to use as high an intensity level as possible. 

MODE 



The mode register is a 3-bit register whose contents determine the way in 
which the next data word will be interpreted. The eight different data 
word foimats (modes) will be discussed in the following paragraphs of this 
section. 



DATA WORD FORMATS 



PARAMETER MODE (000 2 ) 
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Parameter mode is the control mode of the display. A parameter word can 
be used to change the mode, scale, intensity, light pen and/or interrupt 
parameters of the display 



POINT MODE (001 2 ) 
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Point mode is used to specify an X or a Y coordinate location on the dis- > 
play. It can change the mode, light pen, and intensify parameters. A 
data word is interpreted by the display as point mode if the mode register 
was set equal to 001 2 by the previously interpreted word. 



RASTER MODE (010 2 ) 
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CHARACTER MODE (011 2 ) 



1st Character 
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In character mode the display interprets each word as containing three 
alphanumeric characters. Each character is specified by a 6-bit modi- 
fied ASCII code. The display remains in the character mode until an 
escape code is encountered, then the display returns to parameter mode. 
A data word is interpreted by the display as character mode if the mode 
register was set equal to 011 2 by the previously interpreted word of if 
the previous data word was in character mode and did not contain the 
escape character. Bits 0-5 are interpreted as the second chara:ter and 
bits 12-17 are interpreted as the third character. 



VECTOR MODE (100 2 ) 
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In vector mode the display interprets each word as containing vector size 
and direction, intensify, and escape information. The display remains 
in vector mode until the escape bit is set, at which time it returns to 
parameter mode. If the edge of the raster is violated, a flag is set which 
causes a computer interrupt. A data word is interpreted by the display 
as vector mode if the mode register was set equal to 100 2 by the previously 
interpreted word or if the previous word was in vector mode with the 
escape bit equal to 0. 



VECTOR CONTINUE MODE (101 2 ) 
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In vector continue mode the display interprets the word as containing vector 
directions, and intensify information. The vector is drawn from the 
starting point to the edge of the raster. When the vector violates the 
edge of the raster, the display returns to parameter mode. A data word is 
interpreted by the display as vector continue mode if the mode register 
was set equal to 101 2 by the previously interpreted word. 
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INCREMENT MODE (llty 



Esc 


Int 


1st Point 


2nd Point 


3rd Point 


4th Point 
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Right 
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Y 1 Down 
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When in increment mode, the display interprets each succeeding word as con- 
taining information to plot four successive spots; each adjacent to the 
preceding one. A spot can be placed into any one of the eight adjacent 
locations at each movement. A data word is interpreted by the display as 
increment mode if the mode register was set equal to 110 ? by the previously- 
interpreted word or if the previous word was in increment mode and the 
escape bit was not equal to 1. The display remains in increment mode until 
the escape bit is set equal to 1 or it moves a spot past the edge of the 
raster. 
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When in this mode, the display interprets the next word as a jump instruction 
to some location in memory. The subroutine word sets the mode of the next 
word to be interpreted and allows the display of data from nonconsecutive 
memory locations. 



SYMBOL CONVENTIONS 



In the subprogram descriptions , floating point variable names have been 
assigned in accordance with the FORTRAN IV conventions. If integer 
variables are required, names beginning with I, J, K, L, M, or N are 
used. The actual floating point number (or integer) may be used in the 
argument lists in place of a floating point (or integer) variable name, 
if the argument represents input to the subprogram. Constants should 
never be substituted for argument names that represent output from the 
subprogram. 



Details about anything in this section can be found in references [6] , 
[5], [3], and [7]. 



^r 

\ 



SECTION III 



CONTROL ROUTINES 



S ETTING UP ERROR POINTERS: SETSEV 

At the outset of every job, this routine should be called to set up system 
error pointers. 



CALL SETSEV 



CHAIN FILE COMPATABILITY 

The OLMGS has a named common block in which it communicates with itself and 
other systems. This common block is designed for use by CHAIN. A description 
of the common block is in figure 3-1. 



INTENSITY SELECTION: BRITEV,FAINIV,ITSCy 

GRID1V ensures that the bright intensity mode is on. Normally, this intensity 
mode should be left on, since experience has shown that it produces the best 
results. If the programmer wants to change this setting to the faint mode, he 
can use the following statement: 

CALL FAINTV 

The faint mode will always be a setting of 4. Then, to restore the bright 
intensity mode, he can use the statement: 

CALL BRITEV 

Bright mode will always be a setting of 7. To set the intensity and scale to 
any value use the following statement: 



CALL ITSCV Cl, J) 



where I = intensity 0 < I < 7 
J = scale 1,2, 4, 8 






o 

ERIC 
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COMMON /DLT ST/LIST , ZDDPTA , DXXSYY , XXXXDD , YYYYDD , ACDD , AYDD , BXDD , BYDD , DHI GHD , 
DHIGH1 , DLOWDD , DLOWD1 ,MLDD ,MRDD ,MBDD ,MTDD , WIDEDD ,HIGHDD , DDDDYY , DXXDYY , 
DDDSDD , HOLDDD , CAMVDD , YTOPDD , YREGDD , DO , D1 , FRMCNT , SCFL , NOCOMP 
COMMON/ INSCAL/ INTENS , ISCALE , ICHASZ , LENGTH 

FIGURE 3-1 



SCALE SIZE CONTROL: BIGV, SMALL V, CHARV 

The scale may be selected by the programmer. The program statement. 
CALL BIGV 

will set the scale to its maximum setting; and the statement, 

CALL SMALL V 



reduced the scale to its median setting. Character size can be set by calling 
CHARV. The calling sequence is: 

CALL CHARV (k) 

where k is the character size. 



> 

•■Jb 



A 
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SECTION IV 






GRAPHS 



GRAPHING DATA: KWKPLT 



The purpose o£ this routine is to provide the programmer with a quick 
look at the relationship between two variables. KWKPLT will automati- 
cally provide the programmer with a series of linearly connected 
points on a scaled linear grid with identification printing. 

It is not necessary to arrange the coordinates in an increasing or 
decreasing order of magnitude. If the table of X-coordinates are not 
in ascending order, KWKPLT will rearrange them in ascending order 
within the table. The Y-coordinates will be arranged accordingly. 

The calling sequence with identification printing is: 

CALL KWKPLT (X,Y,N,LH,LV, ID) 

where 

X = starting location of a forward stored array of 
floating point numbers representing the X- 
coordinates . 

Y = starting location of a forward stored array of 
floating point numbers representing the Y- 
coordinates . 

N = number of points to be plotted. 

LH = 18 character identification for the X-coordinates. 

LV = 18 character identification for the Y-coordinates. 

ID = 1 label axis. 

= 0 no label on axis. 

The printing routine assumes a full 18 characters 
including blanks. 

PDP340 

This routine provides the programmer with a quick look at the relation- 
ship between two variables. PDP340 will automatically provide the 
programmer with a series of linearly or non-linearly connected points 
on a scaled linear- linear, log-linear, linear-log, or log-log grid with 
identification printing. The calling sequence is: 
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CALL PDP340 (NPLOT , MODE , NCHAR , NPTS ,X,Y ,XMIN,XMAX , YMIN , YMAX ,XLABEL , 
YLABEL ,GLABEL , IERR) 



where 

NPLOT 

MODE 



NQIAR 

NPTS 

X 

Y 

XMIN- 

XMAX 

YMIN- 

YMAX 

XLABEL 

YLABEL 

GLABEL 

IERR 



= 1 eject and start new graph 

= 2 use same graph 

= 1 Linear -Linear 

= 2 Log-Linear 

= 3 Linear -Log 

= 4 Log -Log 

= Character Number 

= Character Number 

= NPTS of X-coordinates 

= NPTS of Y-coordinates 

= Minimum and Maximum of X-coordinates 

= Minimum and Maximum of Y-coordinates 



= 72 character X axis label 

* 72 character Y axis label 

= 72 character heading 

= 1 Normal 

= 2 or 3 Error 
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SECTION V 



> 

GRIDS 



GENERATING A GRID: GRID1V 



In many ways, plotting on the OLMGS is very much like plotting on a sheet 
of graph paper, but there are also distinct differences. For one thing, 
the programmer must create the grid; the picture is completely blank to 
start with. 

Although every line of the grid must be specified on the OLMGS, there are 
advantages to this situation. A hand -plotted graph must be adapted to 
some preprinted foim; more frequently than not, this means that some plot- 
ting area must be sacrificed in order to use the most convenient scale. 

On the OLMGS, the programmer can select a scale that will be easy to read 
and that will accommodate the entire range of data. The number of light 
grid lines, the number of emphasized grid lines, and the spacing between 
lines can be chosen to suit the plot. The prograironer is not restricted 
to the use of a single form for a variety of plots. For each graph, a new 
grid can be tailored to the data. 

The easiest way to create a grid is to call the GRID1V subprogram. At the 
outset, GRID1V makes certain that the Bright Intensity Mode is on. 

GRID1V will produce a grid which has some lines emphasized and some lines 
labeled. Margin space (which may be used for titles) will be reserved at 
the top, left side, and bottom of the grid. Normally, the title margin 
spaces are 2^ raster counts wide. 

Upon completion of GRID1V, scale factors will have been established and 
made available (internally) for the conversion requirements of other 
subprograms; i.e., the conversion of floating point coordinates into 
raster coordinates. 

The call statement for GRID1V appears below, with a description of the 
arguments . 



CALL GRID1V (L, XL, XR, YB, YT, DX, DY, +N, +M, +1, 



+ J , 



1 
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+NX, +NY) 



This integer argument controls the label margins 
computations . 

L f 1 No label margins will be computed. 

L = 1 Compute label margins. 

Floating point values of X for the left-most and 
right -most limits of the grid. 

Floating point values for the bottom 1 imit and the 
top limit of the grid. 

After margin space for titles and labels has been 
reserved, the limits of the remaining space are 
assigned the data values given for XL, XR, YB, and 
YT. Scale factors are then computed; they will 
remain in effect until another GRID1V statement 
is made (or until other action is taken to compute 
new scale factors) . 

Floating point data increment at which vertical 
(specified by DX) and horizontal (specified by 
DY) grid lines will be displayed. If 0.0, no 
lines will be shown. 

Positions are stepped off in DX increments in the 
positive and negative directions from X = 0, and 
in DY increments in the positive and negative 
directions from Y = 0. 

Fixed point integers that cause every Nth vertical 
grid line and every Mth horizontal grid line to 
be retraced for emphasis. If N (or M) is zero, no 
vertical (or horizontal) lines will be emphasized. 

To force the grid to be square, a negative sign 
should be used on N and/or M. (If either N or 
M is zero, the negative sign should go on both 
N and M.) 



I, J Fixed point integers which cause every Ith vertical 

line and every Jth horizontal line to be labeled. If 
I (or J) is zero, no vertical (or horizontal) lines 
will be labeled. 

If I and J are positive, the line labels will lie along 
the X = 0 and Y = 0 lines, provided these lines are 
within the grid limits. If X = 0 (or Y = 0) does not 
fall within the grid limits, labels will be placed in 
a space reserved at the left (or at the bottom) of the 
grid. 

Negative signs can be used on I and/or J to force 
labels outside the grid. Label space is reserved 
at the left if I is negative, or at the bottom if J is 
negative, and labels will be placed in these reserved 
spaces. Note that label margin space is in addition 
to the margin reserved for titles. 

NX, NY Fixed point integers indicating the number of char- 
acters to be displayed in the labels of vertical and 
horizontal lines. 

+NX, +NY The labels will be in a decimal format 
similar to the F-type format. In speci- 
fying +NX and +NY, a decimal point must 
be counted as one of the NX or NY char- 
acters, but the sign is not counted. The 
largest number of digits permitted is 6 
(or 7 if one character is a decimal point). 

-NX, -NY The labels will be in^cientific notation. 

(Example: 1.25x10 .) NX indicates 

the number of significant figures in the 
labels of vertical grid lines, and NY 
indicates the same for the labels of hori- 
zontal lines. The sign, decimal point, 
and exponent will be displayed in addition to 
NX (or NY) characters. NX (or NY) must 
not be greater than 6. 

Examples of GRID1V Usage 



Figures 5-1 through 5-9 are examples of the effect of the various parameters in the 



GRID1V call statement. These examples are reproduced from DEC 340 output. The 
call statement to produce each graph is printed in the figure with the graph. 

Figure 5-1 is a simple grid with the x=0, y=0 lines crossing in the middle of the grid. 

The numeric labels have been placed along the x=0, y=0 lines. For simplicity in the 
illustration, constants were used in the parameter list. In actual usage, variable 
names may be substituted for any parameter. 

Figure 5-2 is similar to Figure 5-1 except that the XL, XR and YB, YT have been 
reversed to show that the scaling routines have no difficulty handling data v/hich 
decreases from left to right and bottom to top. 

Figure 5-3 illustrates the effect of negative values for I, J in the parameter list. 

Note that the numeric labels are outside the grid and that the margins have been in- 
creased to accommodate the labels. 

Figures 5-4 and 5-5 show the same grid with labels in integer notation and scientific 
notation. 

Figure 5-6 has been double exposed to show the effect of negative arguments at N or 
M. The outer frame was produced by the first call statement with the positive argu- 
ment. for N and M. The grid utilizes the maximum available space in both directions 
and is taller than it is wide. The second call statement with negative N and M forced 
the frame to be shorter in the Y direction in order to be square. It is important to > 

have a square grid when representing geometric figures such as a circle or a square. •> 
Note that the first parameter of the second call statement is a 2 which inhibits the 
margin calculation. 

Figure 5-7 illustrates the use of the routine DXDYV to compute some of the values 
for the GRID1V parameter statement. DXDYV is explained on page 5-13. 

Figures 5-8 and 5-9 show the influence upon the grid of the density factor used by 
DXDYV. For the case of Figure 5-8, a density factor of 8. 0 was specified as the 
8th argument of DXDYV. A larger factor, 20. 0, caused DXDYV to derive values of 
DX and DY such that the grid in Figure 5-9 is less dense. 
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Figure 5-1 




5-5 





Figure 5-2 
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400 .00 




CALL MIDtV 0,-400 .0 , 400 .0 ,-400 .0 , 400 .0 , *5.0, *5 .0 , 4, 4 , 9, • ) 



Figure 5-3 
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■ 1000 MOO WOO 4000 1000 0000 7000 

CALL 40 10 1 ▼ ( I i 0.0 , 7000 . 0 , 0 .0, 70000 .0, 100.0 , 1 000 .0 , 5 , 5 , 10 , 10 , 4 , 5 ) 



Figure 5-4 
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I .0X10 



4 .0X10 



I .0X10 



CALL tj* 1 0 1 V < If 0.0 , TO 00. Of 0.0, TOO 0 0.0, 100.0, 1000.0, 5, 5, 10, 10, ~2 » “2 ) 



Figure 5-5 



*«rv 

<> 



O 





too .00 
too .00 



(1 I 50.0 , 1 00 .0 , 50.0 , 
( I , 50 . 0 , 100 . 0 , 50 . 0 , 



100 . 0 , 
100 .0 , 



, 1 , 2 , 2 , 5 , • 

• -1 • 2 . 2 . 3 , « 



too 



Figure 5-6 
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call tmoiv 
call m 101 v 



30 .0 

90.0 



50 . 0 , 1 
50 .0 ,-! 



} 

) 



• 0 




CALL D KD T V (2, 
CALL UUD1 V ( 1 



0.0, 10.0. Dr, M, J , N 
I 0.0, 20 .0 , 0.0, 90.0, 



r i io.o , 
ox, or, 



IERR2 ) 

N , M , ! , 






Figure 5-7 






CALL 0 XD Y V { 1, 0,0, -303.0, Ox, N, ], NX, 9.0, 1ERR J 
CALL OXOYV ( 2 i 93.3, 14.0, Dr, M, J , NY, 9,0, 1ERR > 

CALL WU01V < 1, 0,0, -303,0, 93,3, 14,0, OX, Dr, N , N, 1, J, NX, NY ) 



Figure 5-8 




CALL oxorv ( 1 , 0 . 0 , - 503 . 0 , Ox, N, I, NX, 20 . 0 , I EUR > 

CALL oxorv ( 2 , 93 . 3 , 14 . 0 , OY, M, J, NY, 20 . 0 , IERR ) 

CALL tRIDIV ( 2 , 0 , 0 , - 303 . 0 , 93.3 , 14 . 0 , OX, OY, N, M , I, J, NX, NY ) 



Figure 5-9 
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COMPUTATION OF GRID IV ARGUMENTS: DXDYV 



It frequently happens that the programmer does not have sufficient advance informa- 
tion about the range of data his program will encounter to be able to assign practical 
values to all arguments of GRID1V. In this case, a series of FORTRAN statements 
can be used to determine the upper and lower X and Y bounds. For example, the 
values of XL and XR for a block of data, X, can be computed as follows: 

XL =X(1) 

XR =X(1) 

DO 10 J =2, NPTS where NPTS is the number 

of points in the X block of 
data 

XL = MIN1F(XL, X(J)) 

10 XR = MAX1F(XR, X(J)) 

A similar group of statements can be used to compute YB and YT for the Y block of 
data. 

Once XL and XR (or YB and YT) are known, the routine DXDYV is available to com- 
pute arguments for line spacing, line emphasizing, and line labeling. Two call state- 
ments are available, one for the X direction and one for the Y direction. They are: 

CALL DXDYV(1, XL, XR, DX, N, I, NX, DC, I ERR) 

CALL DXDYV (2, YB, YT, DY, M, J, NY, DC, IERR) 

On each entry to DXDYV, four arguments are furnished by the programmer: 

The first argument is a 1 or a 2 , to indicate whether DXDYV is 
being applied in the X direction or in the Y direction. 

XL and XR (or YB and YT) are defined as in the summary 
of GRID IV arguments. 

DC represents a floating point quantity which limits the density of the 
grid. The grid lines drawn by GRID1V, using arguments furnished by 
DXDYV, will be no closer than DC raster positions. DC should never 
have a value less than 3. 0; values of 8. 0 to 20. 0 are recommended. 

The remainder of the arguments are variables to which DXDYV will assign values. 
Any value previously assigned these variables will be destroyed during execution 
of the subroutine. NEVER USE CONSTANTS FOR THESE ARGUMENTS. 



IERR is an error indicator. It is set to zero if a reasonable grid can be drawn, and 
to one if the parameters given would result in an impossible grid. After execution 
of DXDYV, IERR should always be tested before proceeding to draw the graph. 



In using DXDYV, it should be noted that no provision has been made for generating 
labels in scientific notation. If this is desired, it is necessary to assure that there 
is sufficient space for the longer labels and also to change the sign of NX and NY to 
be negative. 



GRID IV CONTROLS 



Certain features of the basic linear GRID1V can be altered by subprograms that con- 
trol its internal operation. The subprograms can be classified as "set" and "retrieve" 
routines since they permit information to be set by the programmer and retrieved 
during execution of GRID1V. 



The routines that furnish values different from those normally employed by GRIDlV 
are: 



SETMIV, which allows the programmer to make nonstandard 
margin assignments. The companion routine called 
by GRIDlV to retrieve margin values is SETMOV. 

SETCIV, which makes it possible to provide extra space for 

grid line labels. The companion routine is SETCOV. 



Routines that furnish indicators recognized by GRIDlV as signals to execute alternate 
branches are: 



HOLDIV, which assists in holding margins from graph to graph. 
HOLDOV is called to retrieve the indicators. 

SMXYV, which enables the programmer to select a non-linear 
mode of operation. The companion routine is MSXYV. 
These two routines are described under "Log and Semi- 
log Plotting. " 
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Grid Margin Variation: SETMIV, SETMOV 



As discussed in an earlier section, GRID1V normally reserves a strip, 24 raster 
counts in width, at the top, left, and bottom of the grid, for the display of titles. For 
the many applications which require special margin widths, the subprogram SETMIV 
can be called to change the basic specifications. 

One obvious application of SETMIV is to provide margin space for multiple lines of 
printed titles and headings. In addition, and perhaps even more important, SETMIV 
makes it possible to display more than one graph on a frame, or to display a graph 
with its accompanying text. 

The standard GRID1V margin specifications can be altered by the statement: 

CALL SETMIV (MTL, MTR, MTB, MTT) 

Each argument is an integer which specifies, m raster counts, the 
width of one area to be reserved for a margin. 

MTL Width of area for left margin. 

MTR Width of area for right margin. 

MTB Width of area for bottom margin. 

MTT Width of area for top margin. 

GRID.1V does not necessarily use these exact values for the upper and lower limits of 
X and Y. It guarantees that the reserved space will not be overlapped, assigning 
additional space if required for label margins. After the total margin space has been 
reserved, the remaining area will be used for the grid. 

If SETMIV is never called, GRID1V will use the values 24, 0, 24, 24 as MTL, MTR, 
MTB, and MTT, respectively. To return to a standard grid after the margins have 
been altered, restore the standard margin values by 

CALL SETMIV (24, 0, 24, 24) 

The current values of MTL, MTR, MTB, and MTT can be retrieved by using the 
statement 



CALL SETMOV (MTLL MTRL, MTBL, MTTL) 

where the arguments are variables (never constants) to which SETMOV is to assign 
the current margin values. SETMOV was designed for use by GRID1V to retrieve 
current margin values; the programmer will rarely have reason to call it. 



t 

\ 
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Examples . Figure 5-10 shows three grids with the SETMIV and GRID1V call state- 
ments used to produce them. The grid at the bottom was the first one displayed; a 
_l_was used as the first argument of the first GRID1V statement executed, in order 
to change the film frame. The other two GRID1V statements include a 2_ as the first 
argument, to inhibit margin calculation. 

Note particularly the variation in the raster locations assigned by GRID1V to XL in 
each of the grids. This effect is caused primarily by the differences in the specifi- 
cation of NY (the last argument), which gives the number of characters to be dis- 
played in the labels of horizontal lines. In each case, NY has been assigned a value 
just large enough to satisfy the needs of the grid. For the bottom grid, NY = 1; 
for the middle grid, NY = 3; and for the top grid, NY = 5. Since margin space was 
reserved for labels of different lengths, the positions of the left limits, and of the 
corresponding values of X, vary noticeably. Such a nonalignment is often of no im- 
portance, but if it does matter, the programmer may have to make special provisions 
to force alignment. 

PROVIDING FOR SPECIAL LABEL CHARACTERS: SETCIV, SETCOV 

GRID1V computes the starting location of each label, taking into consideration the 
size of the characters used. If the labels are to be placed outside the grid, GRID1V 
assigns space for them, again taking the character size into consideration. Nor- 
mally, labeling is done in CHARACTRON characters (via LABLV). If the pro- 
grammer substitutes a non-system labeling routine for LABLV, it may be necessary 
to furnish adjusted character dimensions to GREDIV. 

To state the dimensions of nonstandard label characters, use 

CALL SETCIV (IW, IH) 

IW An integer which specifies, in raster counts, 
the allowance needed for the width of each 
label character. 

IH An integer which specifies, in raster counts, 

the allowance for the height of a label character. 

If SETCIV is never called, the indicator table contains IW = 8 and IH = 10. 

Obviously, if it is called, the arguments must be compatible with the size of the 
characters employed by the LABLV subprogram used. 

GRID1V retrieves the values of the indicators by using: 
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CALL SETCOV (IWL, IHL) 



NAItttN VARIATIONS FOR tRtDlV 



* 




.i. I 1 L - * » J J 1 i I I J I 1 1 1111 1 i 1 I I 1 I 1 J 

• to to SO 40 SO 

Call ttTMtv ( > 4 , o, Tit, f»> 

Call tRtOiv <t, o ,o , so.o, ioooo.o, aoooo.o, i.o, 2000., 5 , s, 10, 3 , 2, s> 
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CALL SCTNIV ( 
CALL RRtDlV < 



0 * 
0,0 



10 

34*. 3T1) 

, 30 . 0 , 100 , 0 , 



400 , 0 , 1 , 0 , 20 , 0 , 



10 , 3 , 2 , 



-u 



I i I 

t 10 

call sc tni v i t* , o , u , ▼ 1 9 > 

CALL SRI 01 V ( 1 , 0 . 0 , 30 . 0 , 1 . 0 , 



4 . 0 , 1 . 0 , 0 ,t , 



10 , 5 , l, 1 ) 



Figure 5-10 
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The width will be retrieved from the table and stored in the fixed point variable loca- 
tion IWL, and the height will be similarly stored in IHL. (The arguments must not 
be constants.) Note that GRID IV uses this information to control the space that will 
be reserved for labels; it does not control the size of the label characters themselves 
in any way. 

HOLDING MARGINS FROM GRAPH TO GRAPH: HOLDIV, HOLDOV 



For a large graph, it may be necessary for the programmer to display segments of 
the graph in separate frames, and join the segments "tile fashion" to form the com- 
plete plot. If the graphs are to have the same scale, certain equalities should exist. 

1. The range of X in each segment should be equal, and the range 
of Y should be equal. In other words, the quantity (XR - XL) 
should be the same in each segment, and, similarly, the 
quantity (YT - YB) . 

2. The dimensions of the scaled area should be the same from 
segment to segment; that is, the dimensions of the space 
between margins should be equal. 

The programmer can easily provide for equality in the ranges of X and Y, but he 
cannot so readily ensure equality in the scaled areas. Since GRID1V computes label 
margins (and, therefore, total margins) to suit the needs of each graph, the dimen- 
sions of the scaled area may vary. 

One method that will usually give equality of scaled areas is to specify the option 
that forces labels to be placed outside the grid, and to always request the same 
number of label characters (NX, NY) for each segment. If this is not practical, a 
"holding" feature is provided. 

GRID1V can be instructed to hold the label margin spaces used for the preceding 
grid and use them in computing total margins for the next grid. The statement to 
be used is 



CALL HOLDIV (NH) 

If NH / 0, the label margins from the preceding grid 
will be used again. If NH = 0 (as is the case if HOLDIV 
is never called) , label margins will be computed in the 
normal manner. 
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The status of this indicator is tested in GRDD1V by using 
CALL HOLDOV (NHL) 

The value of the indicator will be retrieved and stored in 
the location named in the argument. 

The "hold" may be released by executing the statement 

CALL HOLDIV (0) 

Figure 5-11 is similar to Figure 5-10 except that NY = 5 on all three grids, permitting 
the left limits to be in line. If XL, NX, and NY are equal from grid to grid, the de- 
sired alignment will usually be achieved. 

Figure 5-12 shows four graphs on a single frame. The SETMIV statements used to 
produce the margins for each grid are shown. The programmer must remember to 
set the first argument of the GRID1V statement to 2 so as not to calculate margins . 

Figures 5-13, 5-14, and 5-15 show additional examples of special effects which can 
be obtained with GRID1V when the routine HOLDIV is used to retain grid margins 
from one grid to another. The labeling is self explanatory. 

OPERATIONAL DETAILS OF GRID1V 



GRID1V is, in many respects, an executive routine. It examines the information 
furnished by the argument list and by certain external subprograms, makes decisions 
based on this information, and then calls other subprograms to calculate margins, 
compute scale factors, generate the grid, etc. 

Initially, GRID1V uses BRITEV to ensure that the bright intensity mode is 
"on.'' 



GRID1V then checks certain internal locations to obtain basic information, by calling 
the following subprograms: 

SETMOV Retrieves margin assignments. The programmer 
may have changed the standard specifications by 
a CALL SETMIV. 

SETCOV Retrieves character size specifications that programmer 
may have changed by a CALL SETCIV. 
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CALL sctmiv 
CALL tit 10 IV 



10 

( k'«. 0 , Tit. 22) 

(C| 0.0, 50.0, 10000.0, 



40000.0, 1.0, 2000., 3, 3, 10, 3, t, 3 > 



400 ,0 




CALL tRIDlV ( t, 0.0, 30.0, 100.0, 400.0, 1 .0, 20.0, 5, 5, 10, 5, 2, 5) 




• io ca io 40 tp 

CALL 4CTM1V ( 14 , , o, 24, ?13> 

CALL MIOlV ( 1, 0.0, 90.0, 1.0, 4.0, 1.0, O.2., 9, 9, 10, 3, 2, 9) 



Figure. 5-1 1 
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CALL It TM t V < 24 , 9 34 ( 9 34 , 24 ) 







CALL UTMtV ( 5 3*i 2 2 * 9 31 , 24 ) 








Call SCTMIV < 24 , 9 34 , 24 , 9 3 « > CALL SCTMIV ( 3 3 1 , 22 , 24 , 3 3 « ) 



Figure 5-12 




5-21 

1 1 3 



special tRtCIV EXAMPLES 

HOLD t X AS AN A I*D tN MAINTAINING THE SAME SCALE FROM OME SRIO TO'ANOTHER 



T . 




CALL GRI01V < 1. -10.0. 70.0, 0.0, 7.0, 1.0, 0.1, 3, 3, 10, 10, Z, Z ) 

CALL HOLOIV(l) 

THE HOLO STATEMENT WILL CAUSE LABEL MARGINS TO BE HELD FOR SUBSEQUENT CHIOS 
UNTIL CALL HOLOIVCO) RELEASES THE HOLO 



Figure 5-13 
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Call tMOlV ( 1 ., 70.0, 130.0, 0.0, 7.0, 1.0, o.l, 3, 5, 10, 10, 3, 2 > 

NOTC...NOT ONLY MAS HOLOI* BEEN USED BUT ALSO XL»XR|TB,TT HAVE BEEN SPE C t F I EO SUCH THAT 
*R-*L aNO TT-TB OF ONE 6RI0 EQUAL THEIR COUNTERPART IN THE OTHER 6*10 



Figure 5-14 
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SRECML BRI01V C*AhH.E s 



U|t or HOLOJV WITHIN A 79A MC 



CALL ItTMIV ( 90, 373 i 

cm brioiv ( i* >10.0, 90.0, o.o, 9.0, 9.0, 1.0, «, 2, >4, -2 , 2, 1 » 

CALL HOLD IV (It 

CALL tCTMIV ( 3(4, 91, 30, 373 I 

CALL BRIOIV ( t ( 90.0, 190,0, 0.0, 9.0, ) ( 0, 1.0, 4, 2, -4, >2, 3,1 I 

NOTE... NOT 044 L T MUST XR-XL AND TB-TT X ONE BR I D EtUAl THE CORRE 1*0440 I N 6 BU AN T I TT IN THE 
OTHER BRIO* BUT ALSO THt ARBUMENTt OT IE TN 1 V MU|T HAVE A SPECIFIC 

RELATIONSHIP... |AT THE ARBUMtNT UtT II CM TL ,M TR . M TB ,N TT) . THEM THE 

VALUES 1 023- <MTL*MVR» 'Hi 102 3- IN TM**T> Of ONE BRIO MUST E 9U A L ThE 
CORRE SPOND INC RUANTItT ’.** TK OTHER B4^ 





Figure 5-15 
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MSXYV Determines linear or nonlinear mode. Indicators may 
have been set to nonlinear by a CALL SMXYV. 

HOLDOV Determines whether margin specifications were held 
over from preceding graph. Indicator may have been 
specified by a CALL HOLDIV. 

Computations are made to determine the raster positions to be used as grid boundaries. 
ERRLNV and/ or ERRNLV are called to check the boundaries and the data limits to see 
if a grid can be produced from this information. Finally, XSCALV and YSCALV are 
called to compute scale factors, and LINRV (in the linear mode) and/or NONLNV (in 
the logarithmic mode) are called to generate the grid. 

DETERMINING GRID BOUNDARY POSITIONS 



GRID1V sets aside the margin space specified by a prior call to SETMIV (or the 
standard margin space if SETMIV has not been called). With the exception of the 
small area in the upper right corner used for frame identification, GRID1V does no 
writing in these basi'C margins. 

Then GRID1V tests to see if labels are to be placed (or may extend) outside the grid 
area. Space required for such labels is computed, taking into consideration the type 
of label (fixed point or scientific), the number of label characters specified in GRID1V 
arguments plus space for a sign, and the height and width of the label characters. 

If any space is needed for labels at the left, right, bottom, or top of the grid, it is 
added to the basic margins specified by SETMIV table, to produce the total margins. 
(HOLDIV can alter this procedure by causing label spaces "held" from a previous 
grid to be added when computing the total margins). 

If the option for a square grid is specified, the right, or top, total margin will be ad- 
justed so the remaining area will be square. 

If the total margins and the grid limits (XL, XR, YB, and YT) meet certain error 
tests, they are used as arguments of XSCALV and YSCALV in computing the scale 
factors required for generating the grid and plotting on it. 

Since so many items influence the margin assignments, the grid boundaries will rarely 
fall at the exact raster positions that the programmer might have estimated. If the 
precise raster positions of the boundaries are required in a program, they should be 
derived by converting the limits into raster counts after GRID1V has been called. For 
example: 



IXL = NXV (XL) 



IXR = NXV (XR) 



(The right total margin is 1023 - IXR) 



IYB = NYV (YB) 



IYT = NXV (YT) 



(The top total margin is 1023 - JYT) 



GRID1V ERROR PROCEDURE 

Even if bad input data is used, GRID IV will attempt to produce a grid. The philosophy 
is that some useful information may be revealed, even if the grid is inaccurate. 

The grid limits and total margins are tested by lower-level subprograms, ERRLNV 
(for linear mode) and/or ERRNLV (for log mode). If these tests show that a grid 
cannot be produced from the given information, some data values are manufactured 
so that the program can continue. The manufactured quantities are used only intern- 
ally; data values in the main program will not be affected. (Since ERRLNV and 
ERRNLV are meant to be used only by GRID1V, the call statements are not given.) 

When artificial quantities are used, an error mark (///// ) is placed in the upper 
right corner of the frame by ERMRKV. Although ERMRKV was designed for use by 
GRID1V, the programmer can place this mark on the frame if he uses the statement 



An error mark from GRID1V may indicate that one or more of the following errors 
has been found. 



1. Equal values have been specified for grid limits; that is, 

XL = XR and/or YB = YT. See Figure 5-16. 

2. The specified margins are too wide. In other words, MR + ML 
- 1023, and/or MB + MT - 1023. See Figure 5-17. 

3. In a linear grid, the specified values of DX and/or DY would 
result in grid lines spaced closer than 3 raster counts. 

4. In a log grid the value of one or more of the limits is either zero 
or negative. 

5. In the log mode, there are more than 10 cycles in either the X 
or Y direction. 



CALL ERMRKV 



ERIC i 




I 



tltioiv CftftCM HANDLING 



SPECIAL ttlDtV EXAMPLES 



////////// 






10 



II It 1$ 14 IS 

, 10 . 0 , 10 . 0 , 0 . 0 , 1 . 0 , 0 . 1 , 0 . 1 , 5 , 5 , 10 , 10 , 2 , 2 ) 



it ir io i« 



Figure 5-16 
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CALL tftlOlV ( 1 




MIOIV CftXOM HANOLINt 



»«CUL HIOIV EXAMPLE t 



10 



CALI. itTMIV I «00, fOO t 49^49 ) 

CALL AMIDIV < I, 0.0, 1 0 .0 , 0.0, 10.0, 1.0, 1.0, 3, 3, 10, 10, *, t ) 



Figure 5-17' 




5-28 

12 0 



Examples of GRID IV Characteristics 



Certain special characteristics of the grid generated by GRID1V are shown in the 
following examples. 



Example 1. During construction of a linear grid, vertical lines are displayed at 
intervals "stepped off" from X = 0 in DX increments until the maximum limit of the 
grid is exceeded. Then they are "stepped off" from X = 0 in the negative direction 
until the minimum limit of the grid is passed. Similarly, horizontal lines are 
"stepped off" from Y = 0 in DY increments. Because of this method of construction, 
lines will be generated at the grid limits only if XL and XR are integer multiples of 
DX, and if YB and YT are integer multiples of DY. Figures 5-18 and 5-19 illustrate 
the effect of changing DX and DY. Figures 5-20 and 5-21 illustrate a similar effect 
in the negative range of values. 

Example 2. Vertical lines to be emphasized are "stepped off" from X = 0 by incre- 
ments of N-DX, and vertical lines to be labeled are "stepped off" by increments of 
I-DX. Similarly, horizontal lines to be emphasized are "stepped off" from Y = 0 
by increments of M-DY, and those to be labeled by increments of J-DY. This pro- 
cedure means that labels and/or emphasized lines will not necessarily occur at the 
grid limits. Lines at XL and XR will be emphasized only if XL and XR are integer 
multiples of N-DX, and labeled only if XL and XR are integer multiples of I-DX. 

In the same way, lines at YB and YT are integer multiples of M-DY, and labeled 
only if they are integer multiples of J-DY. 

This method of determining the positions of line labels is responsible for the absence 
of labels on the limit lines in Figures 5-19 and 5-21 . Notice in the example that the 
labels for vertical grid lines are adequate; there is no real need to provide labels on 
the lines at XL and XR. However, the labels for the horizontal lines in this illustra- 
tion are not adequate; they demonstrate another factor to be considered in planning 
for line labels. To be meaningful, labels should appear on at least two lines. In the 
example, the use of J = 5, instead of J = 10, would cause labeling of the lines of 
Y = 25. 0, 30.0, and 35.0, giving a scale that can be read easily. 



O 
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IFCCUL tit ID I V EXAMPLE! 



crrccn of iTtffiNt of r delta i n creme n t s mom ztno 



to 10. CO 

CALL tit 10 IV ( I • 7.3 , 32 .3 i 0.0, 31 .0 , I .0, 2.3 , 3 , 2, 3,4, t, 2 ) 



Figure 5-18 




to 40 

CALL LA ID IV I t, H .0 , 49 .0 , (9 .0 , 99 .0 , I.D, I .0 , 9, 9, ID. ID, t, t ) 



Figure 5-19 
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IFCCIAL tlltDIV EXAMPLE! 



THt NtCATIVt 



KaNU or VALUE! II HANDLED IN A lINILAft WAY 
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-io -is -co -n “io 

CALL WM01V ( 1, -7.5, “32.5, -9 .0 , -Sl.O, 1.0, 2.5, 5, 2, 5, 4, 2, 2 ) 



Figure 5-20 




-10 -4* 



CALL MID1V ( C, *15.0, -45.0, -tS.O, -35.0, 1.0, I .0 , 5, 5, 10, 10 , t, C > 



Figure 5-21 
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Example 3. Under certain conditions, GRID1V purposely omits labels. Figure 5-22 
shows how this can happen. The limits of X used for the grid were -99. 9999 and 
99. 9999, and the limits of Y were 0. 0 and 9. 9999. During scaling, these limits were 
rounded, causing them to be treated as -100. 0, 100. 0 and 0.0, 10. 0. 

This process caused the decimal scales of the grid limits to be larger than the initial 
limits, i. e. , the values of NX and NY specified by the programmer were not com- 
patible with the new limits. To avoid erroneous labeling, label values will not be 
displayed if the rounding process causes the decimal scale to be larger than the 
initial, value. 

The system does not adequately provide for all problems of this type. Such compli- 
cations can be avoided if the programmer considers the rounded values of XL and XR, 
and YB and YT, and when specifying NX and NY. 

Example 4. When X = 0 (and/or Y = 0) lies within the grid limits, and I (and/or J) 
is positive, GRID1V places labels along the X = 0 (or Y = 0) line. If there isn't 
enough space for the label between X = 0 and the left limit of the grid (or between 
Y - 0 and the bottom of the grid), GRID IV will provide space outside the grid area 
for labels, just as if negative values of I (and/or J) had been specified. An illustra- 
tion is shown in Figure 5-23. 

Example 5. The number of label characters, NX (or NY), specified for fixed point 
labels should satisfy the largest and smallest label values to be displayed. Normally, 
the quantity should be the sum of: (1) the decimal scale of the largest value of X 
(or Y), (2) the number of fractional positions required in the smallest value, and (3) 
one more position to provide for the decimal point, if required. The total quantity 
may not exceed 6 (or 7 if the decimal point is included). 

Figure 5-24 illustrates how trouble can occur when NX and/or NY are not specified 
properly. Note that the values of NX = 1 and NY = 1 are adequate for the largest 
values of X and Y, but do not allow for the fractional values required in some of the 
labels. In this example, 3^ should have been used for the values of NX and N Y, in the 
top graph. 

For some applications, it may be advantageous to specify values of NX and/or NY 
less than the decimal scale of the largest label value. In the middle grid of Figure 
5-24, the two low-order positions in the labels have been dropped intentionally, by 
using values of NX and NY that are less than the decimal scales of the largest values 
of X and Y. In this case, the resulting graph has been effectively rescaled. 

This procedure must be used with care. In the bottom grid, Figure 5-24, the 
values of NX and NY are so small that the increment between labeled lines is not 
reflected in the labels that are displayed. 
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IPCCUL tfttD W EXAMPLE 8 




-•0 “BO -40 -10 0 to 40 to 00 

CALL 6MI01V ( It -99.9999 , 99.9999 » 0.0. 9.9999, 4 .0, 1 .0, S, 2, 5, 2, 2, 1 > 



Figure 5-22 



"T 

• to 40 00 

CALL tRIDlV ( t, -0.7 , 100.0 , -3.0 , 100.0 , 2.0 , 10.0 , 3 , 2, 10 , 4 , 3, 3 ) 



Figure 5-23 
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call tuioiv < i. iooo.o, ioooo.o, izoo.o, ioooo.o, 200.0, 400. o, 3, s, s, 3, s, 3 ) 




• * 10 10 II II It It IB 

CALL &RI01V ( ti *00.0, 1300.0, 0.0, 2000.0, 3.0, lDD.D, 5, 5, 10, 3, 2, 1 > 



Figure 5-24 
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LOG AND SEMI-LOG GRIDS 



The earlier examples used only linear scaling and conversion. The modal sub- 
routine SMXYV can be used to alter the scaling mode such that scaling and 
conversions will be made in the logarithmic mode. 

The call statement for establishing the logarithmic mode is: 

CALL SMXYV (MX, MY) 

where MX and MY are scale mode indicators that designate whether the logarithmic 
or linear mode is to be used: 



If 



MX i 0, MY/0 Log in X, log in Y 

MX i 0, MY i 0 Log in X, linear in Y 

MX = 0, MY i 0 Linear in X, log in Y 

MX = 0, MY = 0 Linear in both X and Y 

(to restore linear mode) 



At the beginning of each job, MX and MY are zero, so that linear scaling and 
conversion result if SMXYV is never called. 

If the programmer wants to generate a log or semi-log plot , he must call SMXYV 
to set the logarithmic mode before using any OLMGS subroutines that involve 
scaling or conversion. 

Once the log-log or semi-log mode has been set by SMXYV and the scale factors 
have been established, the function statements NXV and NYV can be used to convert 
data coordinates into raster coordinates, in the same manner as shown for linear 
scaling. 

The contents of the scale mode indicators may be retrieved by using the following 
statement: 

CALL MSXYV (MXL, MYL) 

The indicator for the X scale mode will be 
stored in MXL, and for the Y scale mode in 
MYL. 

GRID1V uses this statement to determine what scale mode has been selected by the 
programmer. 
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RESTRICTIONS ON LOGARITHMIC MODE 



In general, once SMXYV has been called, the programmer can use any of the 
routines to generate a logarithmic display. (One exception: DXDYV should not be 
used to generate arguments for GRID1V in the direction in which logarithmic 
scaling is being used. ) 

However, some of the GRID1V arguments are restricted in the logarithmic mode. 
Following is a list of the 13 arguments, with notations as to the arguments affected 
if the mode is logarithmic in the direction affected by each. 



CALL GRID IV (L, XL, XR, YB, YT, DX, DY, 
±N, ±M, ±1, ±J, NX, NY) 



L Controls margin calculations. 

XL, XR Left and right limits of the grid. May not be negative or zero. 

YB, YT Bottom and top limits of the grid. May not be negative or zero. 

DX, DY Should be set to 1 . 0. If less, only the cycle lines will be 
displayed. If greater, no lines will be drawn. 

N, M Will be ignored; however, an argument must be present. A 
negative sign on N or M will force the grid to the square. 

I, J Will be ignored; however, an argument must be present. 

NX, NY Number of characters to be displayed in the labels. 

Generally, line labels will be placed only at the cycle lines. If, however, the grid 
spans less than one complete cycle, each grid line will be labeled. 

No more than 10 log cycles are permitted in the GRID1V system. 

Examples. Figure 5-25 illustrates a plot that is logarithmic in both the X and Y 
directions. Note the arguments for SMXYV and GRID1V in the coding. 

Figure 5-26 illustrates a semi-log grid, with plotting of the data used in Figure 5-25. 
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CALL 3NXrv (1,1) 

CALL fcHJOlV ( I • l.o, 1 0.0 , 1.0 , 1 000.0 , 1.0 , 1 .0 , 
*OIH73 fLOTTEO BT USfNC CALL P01NTV ( K 



1,2,4) 
1*# A LOOP 



Figure 5-25 
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10000 




call smxvv < o , i ) 

call oxorv ( i, i.o, 20 . o, ox, k. !, nx, iz.o, iek* ) 

CALL tMOlV ( 1. 1.0, IO.D I 1.0 , 10000 .0 , OX, 1.0, N, 1, I, 1, MX, 9 ) 



Figure 5-26 
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BUILDING SPECIAL GRIDS; LINRV, NONLNV 




Two of the lower-level modules employed by GRID1V are useful as building blocks 
for building special grids. LINRV may be used to generate only the vertical portion 
or only the horizontal portion of a linear grid. NONLNV can generate only the 
vertical or only the horizontal portion of a log grid. See examples in Figure 5-27 
and Figure 5-28. 

In addition to generating the vertical and horizontal portions of a grid in separate 
operations, these modules offer other special capabilities: 

1. The programmer can control the lengtti of the grid linos. For 
example, the grid lines can be as short as "time tics." 

2. The programmer can exert greater control over the position of 
line labels by specifying a label reference location. 

3. Selected grid lines can be emphasized and/or labeled in the same 
way as with GRID IV. 

4. DXDYV can be used wRh LINRV in much the same way it is used with 
GRID IV. 

Certain subprograms must be executed prior to the use of LINRV or NONLNV. 

The frame must have been advanced, and XSCALV and/or YSCALV must have 
established scale factors. If a change in scale mode is required, SMXYV must 
have been called. 






ERIC 




The call statement for using LINRV to generate a vertical grid is 



CALL LINRV (1, LYREFR, IYMEN, IYMAX, XL, XR, 
DX, N, I, ±NX, IW) 

For a horizontal grid, the statement is 



CALL LINRV (2, LXREFR, IXMIN, IXMAX, YB, YT, 

DY, M, J, ±NY, IH) 

LYREFR Label reference locations. 

LXREFR For a vertical grid, LYREFR is the Y raster 

coordinate that will be used to position the labels 
of the vertical grid lines. The X raster coordinate 
will vary, and will be computed by the subprogram. 



For a horizontal grid, LXREFR is the X raster 
coordinate that will be used to position the label s 
of the horizontal grid lines. The Y raster coordinates 
will vary, and will be computed by the subprogram. 

The character dimensions, IW and IH (below), should 
be considered in assigning LYREFR and LXREFR. 

For example, LYREFR should be chosen to allow at 
least one character height (IH) below the origin of 
vertical grid lines (IYMEN). LXREFR should be a 
raster position that is at least IW* (NY + 1) raster 
counts to the left of the origin of the horizontal grid 
lines (IXMIN). In both instances, some additional 
space should be added to prevent overlapping. 

If scientific labels are selected (-NX and/or -NY), 
provision must be made for 7 additional label 
characters in the computation of LXREFR. About 
5 raster counts extra must be added to LYREFR 
to allow for the raised exponent. 

NOTE: IW* (NY + 1) allows space for a sign p. 

IYMIN, IYMAX Raster positions that determine the origin and end of 

IXMIN, IXMAX each line. For a vertical grid, these integers are 
raster counts in the Y direction. For a horizontal 
grid, these integers are raster counts in the X 
direction. 
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IW, IH 



The dimension to be allowed for each label character: 
IW for the width, and IH for the height. Since LINRV 
uses LABLV to display the labels, and since the system 
LABLV uses PDP-10 characters, the dimensions 
should be IW = 8 and IH =12 . (If the standard values 
have not been altered, these values can be obtained by 
calling SETCOV). 

The remainder of the arguments have the same definitions given for these terms in 
the discussion of GRID1V. 

The vertical portion of a log grid can be generated by the statement 

CALL NONLNV (1, LYREFE, IYMIN, IYMAX, XL, XR, DX, 

N, I, NX, IW) 

For a horizontal log grid, the statement is 

CALL NONLNV (2, LXREFR, IXMIN, IXMAX, YB, YT, DY, 

M, J, NY, IH) 

The arguments XL, SR, DX, N, I, NX, and YB, YT, DY, M, J, NY have the same 
definitions (and restrictions) given for these terms in the discussion of GRID1V for 
the log mode. The remaining arguments have the definitions given for LINRV. 

Figures 5-27 and 5-28 show examples of portions of grids created by LINRV and 
NONLNV, respectively. 
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examples or use or lin*v 




rt LI *SC<LV 
CALI i 1 MRV 



0.0, 95.0, 50, 5 A A > 

I , 5| 5 , 5 35 , 900 , <5.0 , 95.0, 



10 «0 



?o eo 






2 , a , t , B ) 



• 0 



call yscalv 
call linkv 



15.0, ©5.0, 20 , 575 ) 

2, 520, 350, *.023, 15.0, ©5.0, 2.5, 2, 



2 , 12 ) 



40 



30 



to 



Figure 5-27 
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t H A N P L t ft OF U»C OF NOWLNY 



1 0 




CALL BNXW (1,0) 

CALL XfiCALV ( 10,0, 100.0, SO, S44 > 

Call nonlnv < 1. sis, 353 . too, 10,0, 100.0, 1.0, 



5 , • ) 



10 



.CALL ttfXYV (0,1) 

_CALL YtCALV ( 1.0, 10.0, CO • 3 73 ) 

CALL NOMLNV ( C, 3 ( 0 , 330 . 10(3 . 1.0 , 1 0 . 0 , 1 .0 , 1 , 1 . 2 ,*) 



Figure 5-28 




AXIS LINES: XAXISV, YAXISV 



Axis lines XAXISV, YAXISV enable the programmer to use the axis line feature 
of the DEC 340 to generate horizontal or vertical lines. Horizontal axis lines 
started at a specified raster position will be swept to the right to the specified X 
stop point in raster count. Vertical axis lines started at a specified raster position 
will be swept upwards to the specified Y stop point in raster count. The axis 
generator will not sweep lines down or to the left; therefore, the program auto- 
matically interchanges the coordinates if it is required in order not to stop the DEC 
340 . If no stop point is given in the parameter list, the axis will be swept to the 
right margin or the top mar gin a 

The call statement for sweeping a horizontal line is: 

CALL XAXISV (IX, IY) 

or CALL XAXISV (IX, IY, NSTPT) 

IX, IY The fixed point raster coordinates of the origin of the line. 

IX may have a value from 0 to 959, while IY may have a 
value from 0 to 959. 

NSTPT The fixed point X coordinate of the stop point. 

For a vertical line, the statement is: 

CALL YAXISV (IX, IY) 
or CALL YAXISV (IX, IY, NSTPT) 

IX, IY The fixed point raster coordinates of the origin of the line. 

IX may have a value from 0 to 959, while IY may have a 
value from 0 to 959. 

NSTPT The fixed point Y coordinate of the stop point. 
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Section VI 



SCALING AND CONVERSION 



In earlier sections, scaling and conversion problems have been left to the routine 
GRID1V. However, much of the actual computation is done in lower-level modules. 

This section describes these and some associated modules that provide additional 
tools for some scaling and conversion problems. 

The descriptions may also be useful in clarifying the operation of higher-level routines. 
For example, GRID1V uses XSCALV, YSCALV as a lower- level routine to do scaling. 
Consequently, the comments in this section concerning scaling and conversion equations, 
scale factors, and retrieving and resetting scale factors also apply when the programmer 
uses GRID1V to control scaling. 

Methods of operation in both the linear and nonlinear modes are discussed. The non- 
linear mode built into the system is the logarithmic mode, but the possibility of sub- 
stitution of other nonlinear modes is mentioned. 



BASIC SCALING SUBPROGRAMS: XSCALV, YSCALV 

XSCALV, YSCALV will compute the scale factors for a specified display and store 
them in an internal table for later use by those functions which convert data. The 
calling statements are: 

CALL XSCALV (XL, XR, ML, MR) 

CALL YSCALV (YB, YT, MB, MT) 

XL, XR Floating point values of X for the leftmost and rightmost 
limits of the scaled plotting area. 

ML, MR The amount of margin space to be reserved to the left and 
right of the scaled area, expressed in raster counts (fixed 
point integers). 



-r* 

4*. 



o 
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YB, YT Floating point values of Y for the bottom and top limits of 
the scaled plotting area. 

MB, MT The amount of margin space to be reserved below and above 
the scaled area, expressed in raster counts (fixed point 
integers) . 

XSCALV, YSCALV contains a test for nonlinear mode. If this mode is indicated, XL, 
XR and/or YB, YT will be transformed before the scale factors are computed by the 
basic scaling equations. 



Example 



Figure 6-1 illustrates the relationship of the arguments. The margin specifications 
are: ML = 170, MR = 192, MB = 340, MT = 128. 




Figure 6-1 

XSCALV will assign XL to raster location IX = 170, and XR to raster location IX = 831 
(i.e. , 1023 - 192). YSCALV will assign YB to raster location IY = 340, and YT to 
raster location IY = 895 (i. e. , 1023 - 128). The scaled area will then be the rectangle 
from IX = 170 to IX = 831, and from IY = 340 to IY = 895. 

BASIC SCALING EQUATIONS 



In the following equations, "A" and "B” represent the scale factors computed and 
stored by XSCALV, and "C" and M D" are the factors computed and stored by YSCALV. 
(Since the computation is done in floating point arithmetic, the floating point variable 
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) 



names FML, FMR, FMB, and FMT are used to represent the floating point equivalents 
of the margin values ML, MR, MB, and MT.) 



A _ (1023.- FMR) - FML 
XR - XL 



I 



B = FML - A*XL 



n 



c _ (1023. - FMT) - FMB 
YT - YB 



IB 



D = FMB - C* YB 



IV 



CONVERSION OF DATA: NXV, NYV, IXV, IYV 

Four function subprograms, NXV, NYV, IXV, and IYV, are provided to convert data 
coordinates into raster coordinates. The argument for each of the functions must be 
a floating point quantity; the result will be an integer quantity. 

The following FORTRAN statements show how these functions may be used to convert 
data coordinates X (or Y) into raster coordinates DC (or IY): 

EX = NXV(X) 

IY = NYV(Y) 

IX = IXV(X) 

IY = IYV(Y) 

These four functions are similar in that they all convert data by means of the basic 
equations for data conversion discussed below. They are dissimilar in the way they 
handle off-scale data (that is, data which falls outside the limits XL, XR or YB, YT.) 

The functions NXV and NYV check for off-scale data values. The result IX (or IY) will 
be set to zero if the argument X (or Y) is outside the limits that were used to establish 
the scale. In addition, an error indication is set, as discussed under Off-Scale Error 
Detection. 

The functions IXV and IYV do not test for off-scale data values. The resulting position 
can be outside the plotting area, or even outside the frame, but the value will be properly 
scaled relative to the plotting area. (However, no test is made for the possibility that 



6-3 



130 



the result is greater than 131,071; integer bits above the 17th will be lost.) 

Figure 6-2 illustrates how error testing of the results of NXV, NYV can be used to by- 
pass plotting of points that are off-scale. NXV and NYV were used to convert the points 
along the curve into raster positions, and LINEV was employed to connect the points. 
Since the results of NXV and NYV were tested for zeros, and plotting was by-passed 
whenever a point was off-scale, the curve stopped at the top and right limits of the 
scaled area (outlined). 

In Figure 6-3 two sine curves are shown, one plotted after using NXV, NYV to do the 
conversion, and the other after IXV, IYV were used. As in Figure 6-2, LINEV was 
used to connect the points. NO ERROR TESTS WERE MADE. 

The lower curve shows the line going to zero when off-scale values were encountered 
by NXV, NYV. Note that this curve drops to the bottom of the frame (IY = 0) when 
values that were off-scale in Y were encountered. Also note that the off-scale initial 
and last values of X caused the curve to start and end at the left edge of the frame 
(IX = 0). 

The higher curve was drawn after IXV, IYV were used to convert the points. The curve 
continued past the boundaries of the scaled area when off-scale values were encountered. 
(When IXV, IYV are used, the programmer must decide what action should be taken 
when the result of IXV and/or IYV is < 0 or > 1023.) 

All four functions test for nonlinear mode. If indicated, X (or Y) will be transformed 
before it is converted by one of the basic conversion equations. 

BASIC CONVERSION EQUATIONS 



The following equations show how the conversion functions Convert data coordinates 
X and Y into raster coordinates IX and IY: 

DC = A*X + B 

IY = C*Y + D 

The scale factors A, B, C, and D are those derived from the equations I, II, in, and 
IV. 

Generally speaking, the programmer should use the conversion functions rather than 
writing statements of his own containing these equations. The functions offer the 
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UK or NXV, MTV WITH CMOA ttITINt 



*11 




Figure 6-2 



UK or NXV. NYV i IXV, IYV WITHOUT ERROR TEITINb 




Figure 6-3 
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following advantages: 

a. • They have direct access to the internal table in which the scale factors 

A, B, C, and D are stored. 

& 

b. They check the scale mode, and use a nonlinear conversion if that mode 
is indicated. 

c. NXV, NYV contain a test for off-scale points. 

If it is ever necessary to use these equations directly, the programmer can retrieve 
the scale factors (A, B, C, D) by employing the routine SCLSAV. 

INVERSE CONVERSION UXV, UYV 



The functions UXV, UYV allow the programmer to obtain the coordinates of a specified 
raster location in terms of his data. The following statements show how these functions 
may be used: 

X = UXV (IX) 

Y = UYV (IY) 

Although UXV, UYV represent the inverse of EX = IXV (X) and IY = IYV (Y), the results 
are approximate because truncation occurs in the IXV and IYV functions). 

NOTE 

UXV, UYV CANNOT BE USED 
IN THE NONLINEAR MODE. 

The equations used by UXV, UYV are the inverse of the equations for data conversion: 

x = ix zl b_ 

A 

y _ IY-D 

c 

(The computation is performed in floating point arithmetic. ) 




6-6 

142 



) 




RETRIEVAL OF SCALING INFORMATION: SCLSAV 



The subroutine SCLSAV will retrieve scale factors and other scaling information from 
an internal table and store them in an array named by the programmer. Although 
SCLSAV makes scaling information available for special-purpose conversions, certain 
limit tests, etc. , its principal value is that it permits saving scaling information from 
one program link to another. When a new link is entered, another routine RESCLV, 
can be called to restore the scaling information in the internal table, where it is 
accessible to NXV, NYV, DC/, or IYV. 

The calling statement to retrieve the scaling information is: 

CALL SCLSAV (R) 

R The name of a ten-cell array, named and dimensioned by the 

programmer. 

The storage locations in the block of cells, R, are assigned as follows: 



R(10) 


Minimum IY 


R(9) 


Minimum IX 


R(8) 


Maximum IY 


R(7) 


Maximum IX 


R(G) 


D 


R(5) 


B Scale factors 


R(4) 


C 


R(3) 


A 


R(2) 


Scale mode indicator for Y 


R(l) 


Scale mode indicator for X 
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RESETTING SCALING INFORMATION: RESCLV 



If the scaling information has been stored in COMMON by SCLSAV, it can be reset 
into the internal table when a new chain link is entered. The statement is: 

CALL RESCLV(R) 

where R is the ten-cell array described under SCLSAV. 



NONLINEAR SCALING AND CONVERSION 



The only nonlinear capability built into the system is logarithmic scale mode, used in 
connection with log grids. However, the system design allows the programmer to in- 
corporate some special nonlinear scale mode, by the substitution of one module of his 
own. For this reason, the more inclusive term, "nonlinear, " is used in the following 
discussion, instead of "logarithmic." 



LINEAR^NONLINEAR SCALE MODE INDICATORS: SMXYV, MSXYV 



The subprogram SMXYV allows the programmer to set scale mode indicators. It 
must be called if nonlinear scaling is desired. These indicators are tested within 
other subprograms (GRID1V, IXV, IYV, NXV, NYV, XSCALV, YSCALV). If non- 
linear mode is indicated, additional steps will be taken to handle the selected mode. 
The progr amm er may also retrieve these indicators, using the subprogram MSXYV. 

The statement which sets the scale mode indicators is: 

CALL SMXYV (MX, MY) 



MX ^ 0, MY/ 0 


X nonlinear, 




Y nonlinear 


MX/0, MY A 0 


X nonlinear, 




Y linear 


MX = 0, MY ^ 0 


X linear, Y nonlinear 


MX = 0, MY = 0 


X linear, Y linear 



To reset the indicators for linear -linear, use CALL SMXYV (0,0). If SMXYV is 
never called, the scale mode will be linear in X and in Y, and the internal scale mode 
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indicators will be set as if CALL SMXYV (0, 0) had been executed. 



For chain jobs, note that the values of MX and MY do not carry over from link to link. 
If a scale mode other than linear-linear is desired, it will be necessary to restate the 
SMXYV statement in each chain link. 



The following statement will retrieve the scale mode indicators for testing: 

CALL MSXYV (MXL, MYL) 

MXL, MYL Locations in which the quantities furnished as arguments 
in the last SMXYV statement will be stored . 



NONLINEAR TRANSFORMATION: XMODV, YMODV 



The scaling and conversion equations shown under number conversion apply not only 
to the linear mode, but also to the nonlinear mode if transformed arguments are used. 
The functions XMODV, YMODV are provided to perform such transformation; they 
can be used in statements of the following type: 

XPRIME = XMODV (X) 

YPRIME = YMODV (Y) 

The scaling and conversion subprograms (XSCALV, YSCALV, IXV, IYV, NXV, NYV) 
test the scale mode indicators to see if they have been set to nonlinear mode by a call 
SMXYV. If nonlinear mode is indicated, each of these subprograms will use XMODV 
(or YMODV) to perform a nonlinear transformation on X (orY) before the remainder 
of the scaling or conversion takes place. Since the XMODV, YMODV functions in the 
system compute the log of X or Y, the system nonlinear mode is synonymous with log 
mode. 

For a problem requiring special nonlinear transformation, the programmer can sub- 
stitute subprograms of his own named XMODV and YMODV, with one argument and one 
result. Since the system XMODV and YMODV are physically contained in one sub- 
program, both must be replaced if a substitution is made for either. Obviously, such 
substitute functions must meet the nonlinear scaling and conversion requirements of 
the entire program (or chain link), since they will replace the system functions. 
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System subprograms that use XMODV and YMODV (either directly or indirectly) are: 
GRID1V, KV, IYV, LINRV. NONLNV, NXV, NYV, XSCALV, YSCALV, APLOTV, 
POINTV. 

WARNING: GRID1V WAS DESIGNED ONLY FOR THE LINEAR, LOG, AND 

SEMI-LOG OPTIONS. SUBSTITUTION OF A SPECIAL TRANS- 
FORMATION FUNCTION MAY CAUSE UNEXPECTED DIFFICUL- 
TIES IF USED BY GRID1V. 

OFF-SCALE ERROR DETECTION 



Whenever there is a possibility that off-scale data points might be encountered, error 
tests should be made by the programmer. For example, APLOTV sets an error 
indicator which should be tested. A zero result from NXV or NYV nearly always 
indicates an error; a test should be made for this condition. 

There are additional situations which demand special error detection procedures. For 
one thing, if no left and/or bottom margin space is reserved, the conversion of XL and/ 
or YB can produce a legitimate zero result from NXV, NYV. More important, the 
programmer may be using NXV or NYV indirectly, via other modules, and thus be un- 
able to test the results. For these reasons, additional subprograms are provided for 
detailed analysis of conversion errors resulting from off-scale points. 

Keep in mind that the special procedures which follow are designed for unusual situations, 
in which normal error testing does not suffice. 



SET CONVERSION ERROR INDICATORS: SCERRV 



Two internal cells are used by NXV and NYV to store indications of successful or un- 
successful data conversion. The subroutine SCERRV allows the programmer to assign 
two cells which NXV, NYV will use in place of the internal error cells. In this way, 
the programmer can name error indicator locations that are accessible to his program. 
The call statement is: 

CALL SCERRV (KX, KY) 

When NXV converts a quantity successfully, it will place a "0" in KX; when unsuccessful, 
it will store a "1" in KX. NYV will use the cell KY in the same way. The cells named 
will be used in each subsequent execution of NXV, NYV (including execution via other 
subprograms) until new cells are named by another call to SCERRV (or the internal 
cells are reset at the beginning of a new link of a chain job). If possible, tests of these 



error cells should be made as soon after execution of NXV, NYV as possible, to avoid 
any possibility that they might be altered by subsequent executions of NXV, NYV. 



The following practical examples show how SCERRV can provide error indicator cells 
in connection with POINTV. Since POINTV uses NXV and NYV only once, the named 
error cells will still contain indications of off-scale errors when control is returned 
from POINTV to the calling program. (The examples assume that the scale factors 
have already been established.) 



Example 1 



CALL SCERRV (KX, KY) 



DO 300 I = 1, N 



CALL POINTV (X, Y, NS) 



IF (KX*KY) 700, 300, 700 Test for non-zero KX and/or KY 



300 CONTINUE 



700 CALL DUMP 



Dump when an off-scale point is 
encountered 



yifT V’ 
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Example 2 



CALL SCERRV (KX, KY) 
Ml = 0 
M2 = 0 

DO 300 I = 1, N 
CALL POINTV (X, Y, NS) 
Ml = Ml + KX 
300 M2 = M2 + KY 



In Example 2, Ml and M2 will contain the total number of points that are off-scale. 
The contents of these locations may be printed out at the end of the job or after all 
points have been plotted on each frame. 



SAVING AND RESETTING ERROR INDICATOR CELLS: SERSAV, SERREV 



As has been pointed out, the locations assigned by SCERRV will be used as error 
indicator cells until new ones are named (or a new link of a chain job is loaded). As 
a result, subsequent executions of NXV, NYV, whether on the same level or on a lower 
or higher one, will alter the contents of the error cells (setting them to _l_for an un- 
sucessful conversion, to 0^ for a successful conversion). 

Two problems can occur. The obvious one is that indications of error may be masked 
by a subsequent execution of NXV, NYV. This can usually be avoided if the programmer 




completes error tests before there is any possibility that the contents of the error cells 
might be changed. 

The less obvious problem is that the error cells named in one level of a program will 
not be accessible for error testing in other levels unless special action is taken. One 
way to solve this problem is to carry the error cell names in the call statements of the 
subprograms, not always a desirable solution. 

The modules SERSAV and SERREV, used in conjunction with SCERRV, offer a con- 
venient means for avoiding these problems in multi-level jobs. The call statements 
are: 



CALL SERSAV (LOCX, LOCY) 

This subprogram saves, in LOCX and LOCY, the locations of the cells that 
are currently being used for off-scale indicators. 

CALL SERREV (LOCX, LOCY) 

This subroutine resets the off-scale error cell locations that were saved by 
SERSAV. 

The arguments, LOCX and LOCY, must be variable names (either fixed or floating 
point) which are not used for any other purpose between the execution of SERSAV and 
of SERREV. 



3 
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Section VII 



PLOTTING 



PLOTTING DATA 



The 64 PDP-10 characters are shown in Figure 7-1 along with the selection 
code for each. 

Reduced to fundamentals, printing or plotting of one PDP-10 character in- 
volves the selection and display of that character at a specified 
position on the raster. (A basic subprogram, PLOW, can be used to 
select and display one character at a time.) The higher-level subprograms, 
however, contain features that make each one suitable for a specialized 
purpose: plotting, printing, or labeling. These specialized features 

of the higher-level routines make them appear to be distinctly different 
from each other. 

POINT PLOTTING SUBPROGRAMS 



Since point plotting usually involves the scaled representation of a 
physical phenomonon, most point plotting subprograms accept physical 
data coordinates for position information. During execution of these 
subprograms , the data coordinates are converted into raster coordinates . 
Scale factors must have been established for the plotting routines to 
use in making these conversions; this requirement can be satisfied by a 
prior entry to GRID1V. 

Scale factors are not normally saved from link to link in a chain job. 
Consequently, plotting should be done within the link in which scale 
factors are computed. 

PLOTTING AN ARRAY: APLOTV 



APLOTV was designed for situations in which a large number of X values 
are stored in one array and and the corresponding Y values are stored 
in another array. It is possible to plot the entire set of data with 
one entry to APLOW. If desired, only a portion of the data can be 
plotted . 
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It is also possible to use different PDP-10 characters as symbols to identify 
curves. The programmer furnishes a table (of one or more PDP-10 charac- 
ters) that will be selected cyclically. 

APLOTV also keeps a tally of the number of off-scale points encountered. 



The calling statement is: 



CALL APLOTV (±N, XARRAY, YARRAY, JX, JY, ±NC, MARKPT, IERR) 

N Controls the number of points to be plotted. N is actually the 

number of points if all the data points in the arrays are to be 
plotted in succession. The value of N may be computed by 
letting K in the following formulas equal the number of points 
to be plotted: 

N = K*JX or N = K*JY, whichever is larger. 



Usually, N should be positive, to indicate that the data arrays 
are in normal FORTRAN order of storage. However, if the 
arrays are in increasing order of storage, N should be negative. 



XARRAY Normally, the names of the arrays of floating point data to be 

YARRAY plotted. Since these arguments must name the locations of 

the data coordinates of the first point to be plotted, subscripts 
may be necessary. Both arrays must be in the same order 
of storage. 

JX, JY Fixed point positive integers giving the increments to be 
added to the subscripts of XARRAY and YARRAY as each 
point is plotted. 



NC The number of characters in the array MRKPT to be used as 

plotting symbols (usually the size of the array MRKPT). 
Positive NC indicates that MRKPT is a normal FORTRAN 
array. Negative NC shows that MRKPT is stored in increasing 
locations in core. 



MRKPT The array that contains the plotting 

character(s) to be used. These will be used cyclically, with 
the first one being used again after the NCth one has been used. 



To set up this information when only one character is to be 
used, the argument can be one of the following: 



i 
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a. A Hollerith argument for one FORTRAN character, as 
"1HO. " 

b • The name of a location containingan "0" , as "MRKPT. ' 



a . An array read in by the A format, one character to a 

location. Such an array is not restricted to integer-type 
names ; it could be called PTMRK, for example, if the 
contents are not used in a computation. 

IERR The name of an error location supplied by the programmer. 

Any point that falls outside the grid drawn by GRID1V will not 
be plotted. Instead, a count of such points will be stored in 
IERR. 

APLOTV plots only one character per point. If a center dot is desired, APLOTV 
may be repeated, using the plotting dot as the symbol. The use of APLOTV is 
illustrated by the Figures 7-2 and 7-3. 

In Figure 7-2, it is assumed that the data are stored in two arrays: X(l), X(2), . . . , 
X(25), and Y(l), Y(2), . . . , Y(25). A prior entry was made to GRLD1V, using the 
statement: 

CALL GRID1V (1, 0.0, 14.0, 0.0, 180.0, 0.2, 5.0, 5, 4, 10, 8, 3, 3) 

Note that the X and Y limits cover the range of the data to be plotted. Also, since 
GRID1V was entered first, the necessary Beale factors have been established for 
APLOTV. /The calling statement of APLOTV for this example is the following: 



DATA II ITOfttD XI ,Xt 



KtJ AND T 1 ,T* 



I 

£ 



[ 1 



£_ 



4— 



TT 



i 

r- 



10 . 



• t. 4. I, 9 . 

CALL MIDI* (I, 0.0, 14.0, 0.0, 140.0, O.t, 5,0, 5, 4, 10, 4, 3, 5) 

CALL A P LOT ¥ (t» , X, Y, 1, S, 1, 34, It** ) O* CALL A*LCTV (25, X, Y, 1, 1, 1, 1HO, !C*R ) 



It. 



Figure 7-2 
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Figure 7-3 
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CALL APLOTV (25, X, Y, 1, 1, 1, XT, IERR) 
or CALL APLOTV (25, X, Y, 1, 1, 1, 1HO, IERR) 



N = 25 The array size, 25 in this case, is equivalent to the number of 
points to be plotted. This can only be true when all points in 
each array are to be plotted. 



X, Y Names of the X and Y arrays (properly dimensioned). 

JX = 1 Since every X versus Y is to be plotted, the fourth and fifth 

JY = 1 arguments are set to 1. 



NC = 1 Only one plotting symbol is to be used; hence the sixth argu- 
ment is set to 1. 



MARKPT In the two calling sequences shown, one specifies the plotting 
symbol by a symbol ; the other uses a Hollerith argument. 



Figure 7-3 shows one way that APLOTV might be used to plot a family of curves. 

The example assumes that the XA and YA arrays are properly dimensioned. The 

data is stored XA(1), XA(2) XA(10); corresponding Y's for the three curves 

are stored in the array YA in the order A(l), A(2), . . . , A(10), B(l), B(2), . . . , 
B(10), C(l), C(2), . . ., C(10). PTMRK is the name of a three-word array which 
was loaded as the BCD equivalents (read in by the A format) of the characters O, X, 
and *. 

In the example, A(I), B(I), and C(I) were plotted versus X(I) each time APLOTV was 
entered. A DO loop was used to proceed to the next value of X, so that a total of 10 
entries were made to APLOTV. The coding was: 

DO 1240 I = 1,10 

1240 CALL APLOTV (30, XA(I), YA(I), 0, 10, 3, PTMRK, IERR) 

Note that the plotting symbols are used cyclically, returning to the first one when 
the array PTMRK is exhausted. If desired, a center plotting dot can be superimposed 
upon the plotting symbol by repeating the entries to APLOTV with the plotting dot 
used for PTMRK. 



PLOTTING INDIVIDUAL DATA POINTS: POINTV 

For each entry to POINTV, one symbol is plotted (with or without a center dot). The 
coordinates may be specified as floating point data, which POINTV will convert into 
raster coordinates. Scale factors must have been established; this can be ac- 
complished by a prior entry to GRID1V. 
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The calling statement is: 



CALL POINTV (X ,Y ,+NS ,ANY) 

WHERE X,Y are coordinates of the point to be plotted, 
stated as floating point data values . 

+NS If NS is minus one, there will be no center 
dot plotted. If NS is positive one, a 
center dot will appear. 

ANY any character to be displayed at the point 
pointed at by POINTV. 

Points outside the scaled area will not be plotted. Until the programmer 
learns how to detect off-scale points, he should be sure that the data 
coordinates will fall within the limit of the scaled area. 



An alternate version of this subprogram allows the programmer to specify 
position information in raster coordinates . This is particularly useful 
when the programmer wants to construct a legend in the margin, showing 
the symbols used and their meaning. The alternate call statement is: 



CALL P0INTV(IX , IY , +NS , ANY) 

Where IX, IY are raster coordinates of the point 
to be plotted; fixed point. 

+NS If NS is a minus two, there will be 
no center dot plotted and the routine 
will accept IX and IY as raster co- 
ordinates. If NS is a positive two 
a center dot will appear and IX and 
IY will be accepted as raster coordi- 
nates . 

ANY Any character to be displayed at the 
point pointed at by POINTV. 



BASIC PLOT-PRINT SUBPROGRAM: PLOTV 



Any of 64 PDP-10 characters can be displayed at a specified raster position 
by using PLOTV, the basic subprogram used as a lower-level modu] ' of other 
routines . 
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However, the programmer may find it useful when other plotting or printing sub- 
programs are not suitable. The call statement is: 



CALL PLOTV ( IX , IY , ANY) 

Where IX, IY are fixed point raster coordinates at which the 
character will be displayed. 

ANY is any desired character. (See Figure 7-1). 
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Figure 7-5 
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LINE GENERATION: LINEV 



II NEV connects two points by a straight line composed of vectors, joined end-to-end. 
The arguments for LINEV, which specify the points to be connected, must be given 
in raster counts. As described above, the programmer may connect two data points 
by a line if he first uses the functions NXV and NYV to convert the data coordinates 
into raster coordinates. (If there is a possibility that the data points being converted 
may be off-scale, the conversion results should be tested for errors before LINEV 
is executed.) The calling statement is: 

CALL LINEV (1X1, IY1, 1X2, IY2) 

1X1, IY1 Raster coordinates of one end point. 

1X2, IY2 Raster coordinates of the other end point. 

Figure 7-6 contains an illustration of the use of LINEV. 

LINE2V is used to draw a line from a fixed point in some direction specified by DX 
and DY. The calling sequence is: 

CALL LINE2V <1X1, IY1, IDX, IDY) 

1X1, IY1 Raster coordinates of starting point. 

IDX, IDY Number of raster points that the line is to be extended in the 
X and Y directions. DX and DY will be handled modulo 64. 

In either LINEV or LINE2V, a floating point data value may be utilized, if scaling 
has been established, by utilizing the function subprograms NXV, NYV as follows: 

CALL LINEV (NXV (XI), NYV(Y1), NXV(X2), NYV(Y2)) 

Figure 7-7 is an example of the use of LINE2V. Each line in the figure is produced 
by incrementing IDX or IDY. 
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CALL XSCALV (-511.0, 512.0, 0, 0) 

CALL YSCALV (-511.0, 512.0, 0, 0) 

Z0 =0.0 
Z1 =4.0 
Z2 =4.0 

CALL POINTV (Z0, Z0, -16) 

DO 5 1=1, 63 
INC = I 
JNC = -INC 

CALL LINE2V (NXV (Zl), NYV (Z0), 0, INC) 
CALL LINE2V (NXV (Z2), NYV (Z0), 0, JNC) 
CALL LINE 2 Y (NXV (Z0), NYV (Zl), JNC, 0) 
CALL LINE 2V (NXV (Z0), NYV (Z2), INC, 0) 
Zl = Zl +3.0 
5 Z2 = Z2 - 3.0 




Figure 7-7 
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Section VII 

TITLING AND LABELING 
TITLING AND LABELING SUBPROGRAMS 



The printing and labeling subprograms enable the programmer to affix titles 
and other identifying information to a picture. Three subprograms of this 
type will be introduced: PRINTV, APRNTV, and LABLV. Other means for 

printing and labeling will be given in the section on printing. 



For many applications, the positions of titles and labels must be independent 
of the scale. Therefore, printing and labeling subprograms accept position 
information in raster coordinates. This contrasts with the plotting routines, 
which includes facilities for the conversion of data into raster counts. 

Titles or labels may be positioned relative to data. The conversion functions 
discussed later can be employed to find raster coordinates from floating 
point location data. 



HORIZONTAL TITLES: PRINTV 



This subprogram allows the programmer to display horizontal titles composed 
of characters. The call statement provides for printing characters read in 
by A- type format. The call statement is: 



CALL PRINTV (N,ASCTEX, IX, IY) 

N The number of characters to be printed. 

ASCTEX An array containing the (A- type format) text to be 
printed . 



-Tv 
-4 Tv 



O 
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IX, IY The raster coordinates for the center of the 
first character. 



VERTICAL TITLES: APRINTV 



This subprogram can be used to display vertical titles composed of 
PDP-10 characters. Each individual character will be upright . The 
call statement provides for printing characters read in by the A 
format or characters stored as a Hollerith argument. Since APRINTV 
prints text or character, the spacing of the characters is controlled 
by arguments specified by the programmer. 

The call statement for APRINTV is: 

CALL APRINTV (INCRX, INCRY, N, ASCTXT, IX, IY) 



INCRX Increments used to space the 

INCRY characters in the X or Y. 

direction, given in raster 
counts. For vertical titles, 
INCRX will be zero and INCRY 
should have a negative value. 
(It is suggested that INCRY 
fall in the range between -12 
and -18 for vertical titles 
in most applications) . 



The remaining arguments are as specified under PRINTV. 
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FIXED POINT LABELS: LABLV 



LABLV was developed for GRID1V to employ in labeling grid lines, but will 
prove helpful when the value of computed quantities must be printed. The 
routine performs one chief task: it displays a floating point number at 

the raster coordinates specified. 



The calling statement of LABLV is: 

CALL LABLV (D, IX, IY, NCHAR, NT, NDMAX) 

D The floating point quantity to be printed. 

IX, IY The raster coordinates which will position 
the first character of the label. Note 
that this first character may be a leading 
blank. If the quantity to be displayed is 
negative, the minus sign will be displayed 
one character space to the left of IX, IY. 

NCHAR Number of characters to be displayed, including 
leading blanks and the decimal point, if any. 

NCHAR is limited to 6 (or 7 if one of the 
characters is the decimal point) . 

NT The number of times each character is to be 

displayed (number of over-strikes) . Norma 11 y 
this should be 1, but 2 or more may be chosen 
if a darker label is desired. 

NDMAX Maximum decimal scale; i.e., maximum number of 
characters to be displayed to the left of the 
decimal point. 

An integer quantity may be displayed by first changing it to floating 
point form and then using LABLV. 

An alternate form of LABLV may be used to display labels in scientific 
notation. The call statement shown is used with the following changes. 

NCHAR Number of significant figures to be displayed. 

NCHAR may be less than or equal to 6. The 
negative sign will result in the use of 
scientific notation. 

NDMAX May be any fixed point quantity. Since the 
right adjustment of these labels will not 
be necessary, the value of NDMAX will be 
ignored . 
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NCHAR will affect the format in the following ways : 

NCHAR FORMAT 

1 YlxlO+YY 

2 Y.YxlO+YY 

3 Y. YYxlO+YY 

4 Y. YYYxlO+YY 

5 Y.YYYYxlO+YY 

6 Y.YYYYYxlO+YY 

Since the space required for these labels will be greater than that required 
for the fixed point format, the programmer should allow NCHAR + 7 character 
spaces in width and 1-1/2 spaces in height as a minimum; it may be necessary 
to allow even more to avoid overlapping other images. 




) 
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S— V 



SECTION IX 



PRINTING 



OLMGS 340 PRINTING: 1X540 

The purpose of this routine is to provide the programmer with a quick 
method of printing on the CRT. 

The calling sequence is: 

CALL TX340 (A , NCHAR , NLINE , NCOL , IER) 

where A = the array of character to be printed 

NCHAR = Number of characters. 

NLINE = Number of lines. 

NCOL = Number of columns 

IER = 1 normal 

2 Error 



« 
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SECTION X 



MISCELLANEOUS ROUTINES 



These routines are classified because they are never used directly by the 
programmer; however, many may be used in other programs. 



UNPACKING CHARACTERS: UNPACK 

This subprogram is used as a low-level subprogram to unpack characters. 
Characters in one vector are unpacked into a second vector one 
character per word (A- type format compatible) . The call statement 
for UNPACK is: 



CALL UNPACK (VEC,UVEC) 
where : 

VEC is a vector of characters (any length) . 

UVEC is a vector of length five , which will 
contain the unpacked characters . 



SETTING SYSTEM FLAGS: FLAGS 

System error indicators may be set by the programmer while scaling data, 
by calling FLAGS. Note it is not necessary for the programmer to scale 
any data in this system; however this facility is available if the 
programmer uses this option. The calling statement is: 

CALL FLAGS (I ,N) 



where : 

I is the. value of the indicator. 
N=0 set X indicator to I„ 

^0 set Y indicator to I. 



DRAWING VECTORS: PLTW5 



This subprogram is used to draw a vector from a fixed point in some 
direction specified by IDELX and IDELY. The calling sequence is: 

CALL PLTW5 (1X1 , IY1 , IDELX , IDELY) 

IX1,IY1 raster coordinates of starting point. 

IDELX, IDELY number of raster points that the line 
is to be extended in the X and Y 
directions. 

10-1 



DISPLAY PREMITIVES 



These routines are described in reference [3] . 

PGEN(IX,IY, INTENSITY) 

VGEN( IX, IY, INTENSITY! 

TEXTP 

PARMS (ISCALE , INTENSITY) 
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SUBROUTINE 



PAGES 



SETCOV 5- 
SETMIV 5- 
SETMOV 5- 
SETSEV 3- 
SMALLV 3- 
SMXYV 5- 
TX340 9- 
UXV 6- 
UYV 6- 
XAXISV 5- 
XMDDV 6- 
XSCALV 6- 
YAXISV 5- 
YMODV 6- 
YSCALV 6- 
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The Division of Computer Research and Technology, NIH, will 
issue on an irregular basis technical documents which we believe 
will be of particular interest to the biomedical community. 

These reports will include detailed descriptions of relevant 
computer programs and instructions in their use (as well as some 
theoretical background) , in hopes that interested scientists will 
be encouraged to gain first-hand experience in applying them. 

In some cases, such reports may serve as foci around which DCRT 
will structure training courses to expand the knowledge and 
experience of NIH staff in applying computer science to problems 
of research and management. Circulation of these reports within 
the biomedical community broadly is, of course, encouraged. 
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A. W. Pratt, M.D. , Director, DCRT 
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SECTION I 



f INTRODUCTION 



The basic philosophy of the OLMS is the minimization of storage used 
during the time when display is INITED, while at the same time keeping 
the system flexibility under the control of the user. 

This part is divided into four sections. The contents of the sections 
are listed below: 



r 



l 




SECTION II : 



SECTION III: 



SECTION IV : 



A discussion of the routines used by the 
interpreter for the command language. 

The routines used to generalize the OLMGS 
under INTER are discussed. Routines 
which input data and graphical functions 
using the RAND TABLET and Function Key 
routines are discussed. 

The system I/O is discussed in this section. 
Reasons for writing I/O are: 

(1) To detect as many errors as possible 
before an EXIT to the monitor. 

(2) A magnetic tape directory system was 
needed . 

(3) The magnetic tape directory system 
which existed, called FILER, had no 
user level error returns that could 
be handled before an EXIT to the 
monitor . 
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DIRECTORIES 



The OLMS has two directories called PAGET and SLIBD. PAGET and XLIBD 
are the system file and library directories respectively. In the present 
configuration PAGET and XLIBD are set up for fifty files in each directory. 
Each file has an associated file description of fifty characters. Each 
directory is a total of 700 words. An example of the directories is given 
in figure 1. 



NAME 


EXT 


PROP 


DESCRIPTION 
50 CHARACTERS 



























PAGET 



NAME 


EXT 


PROP 


DESCRIPTION 
50 CHARACTERS 



























XLIBD 
FIGURE 1 



SYSTEM GENERATION 

A subsystem of this system can be generated with a few minor changes 

in some of the routines associated with the interpreter. For example, a 

user could get along with twenty or less files and no file description. 

This saves 1360 words of storage. Clearly, this is a compromise between 

generality, storage and processing time. Many subsystems can be developed 
to fit the individual user’s needs. 
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INTERPRETER 
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INTERPRETING COMMANDS: INTER 

This routine interprets the command language. The string to be 
interpreted must be in a common block called STRBUF. 

The calling sequence is: 

CALL INTER 

No arguments are present since this routine is defined on a common 
block . 

STORAGE ALLOCATION: GETCOR 

Storage is allocated by the interpreter by calling GETCOR. 

The calling statement is: 

CALL GETCOR (N,IER) 

N = The amount of storage to allocate. 

IER = 1 The amount of storage asked for is not available. 

= 0 The storage is allocated. 

Storage is allocated in 1024-word blocks. Any details about how storage 
is allocated can be found in reference (1) • 

TASK ALLOCATION: INTERP 

After each task is complete in the OLMS this routine should be called 
to return control to the OLMS. This call should be made at the point where 
the user wishes to exit or return from his routine. 




2-1 



The calling statement is: 

CALL INTERP (JCHDEV) 

JCHDEV is the name for the device where the chains are stored. 

All modeling systems running under the OIMS must make this call in the, 
instruction sequence. JCHDEV is a variable in a common block called SYSDDD. 

LIST STRING INTERPRETER : LISTT 

The list commands are interpreted by a sub interpreter of INTER. LISTT 
is used to give flexibility. A part of the interpreter may be a chain file. 
This minimizes storage in case a user wants a system with a resident inter- 
preter. 

The calling statement is: 

CALL LISTT 

This routine assumes that the string is in a common block called STRBUF. 

FILE NAME DELETIONS FROM SYSTEM DIRECTORIES : DELEFL 

File names are deleted from the system directories using DELEFL. The 
director'/ is determined by the state of the ITAMAS switch located in the 
common block called COMSWT. The file name deleted is the name equal to 
XNAME, the extension equal to EXT and the property equal to PROP located in 
the common block called COMFIL. LOOK must be called before DELEFL to set 
the I FOUND pointer. 

The calling statement is: 

CALL DELEFL 

ITAMAS = 0 Before call delete from PAGET. 

Before call delete from XLIBD. 
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ITAMAS = 1 



BLANK VECTOR: BLANKS 



This routine blanks a buffer of five words. This buffer is used to 
send diaracters to a name syntax checker. The buffer or vector is called 
XPNAM and is located in a common block called PNAME. 

The calling sequence is: 

CALL BLANKS 

LISTING MODE: MODE 

The OLMS has eight modes. This routine is defined on MMODE. MMODE 
is the system mode indicator. The mode indicator is located in the common 
block called MODTAB. 

The calling statement is: 



CALL MODE. 

DEVICE CHARACTERISTICS: CHARD 

This routine determines the physical characteristics associated with a 
logical device. The DEVCHR UUO is used in this routine. 

The calling sequence is: 

CALL CHARD (DEVIC , IACL , IACR) 

DEVIC A sixbit device name. 

IACL The content of the left -half of the accumulator. 

IACR The content of the right-half of the accumulator. 

Details about this routine can be found in reference (1) . 

FIND ADDRESS: GETADR 

This routine returns the address of the variable in the first argument. 
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The calling sequence is: 



CALL GETADR (VAR , I RVAR) 

VAR Find the address of this variable. 

IRVAR The address of the variable. 

BLOCK TRANSFER DATA: TRANSF 

This routine will block transfer data from the place pointed at by the 
first argument to the place pointed at by the second argument. 

The calling sequence is: 

CALL TRANSF (IF, IT, N) 

IF A link variable pointing to the place where the data 
will be moved from in the transfer. 

IT A link variable pointing to the place where the data 
will be moved to in the transfer. 

N One less than the number of words to be moved. 

The pointers can be found by using GETADR. Data-may be transfered in 
the opposite direction by exchanging the link variables. 

CONVERT TO SIXBIT: CONVER 

This routine converts ASCII characters to SIXBIT characters. 

The calling sequence is: 

CALL CONVER (XNAME) 

XNAME Contains the ASCII characters when the routine is 
called and SIXBIT when the routine returns. 
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LOOK FOR NAME: LOOK 



This routine looks a' the content of the system directories to determine 
if a file is present. The directory is determined by the state of the ITAMAS 
switch located in the common block called COMSWT. The file name looked for is 
the name equal to XNAME, the extension equal to EXT and the property equal 
to PROP located in the common block called COMFIL. After the search is made 
for the file the switch LOCATE in COMSWT is set to one if the file name is found; 
otherwise, LOCATE is zero. If the file is found a pointer called IFOUND is set 
to point to the location in the directory where the file name is located. 

IFOUND is a pointer located in the common block called CQMPNT. 

The calling sequence is: 

CALL LOOK 

ITAMAS=0 Before call will look at PAGET. 

ITAMAS=1 Before call will look at XLIBD. 

LOCATE=0 After the return indicates the absence of the file name. 

L0CATE=1 After the return indicates the presence of the file name. 

IFOUND Pointer to the file name location. 

PLANTING FILE NAMES: PLANT 

This routine is used to plant file names in the system directories. The 
directory is determined by the state of the ITAMAS switch located in the 
common block called COMSWT. The file name planted is the name XNAME, the 
extension EXT and the property PROP located in the common block called COMFIL. 
This routine calls LOOK and TOP. If an error occurs an error message will be 
written before the return. 
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The calling statement is: 

CALL PLANT 

ITAMAS=0 Before call plant in PAGET. 

ITAMAS=1 Before call plant in XLIBD. 

PLANTING FILE DESCRIPTIONS: PDESCR 

This routine is used to plant file descriptions in the system directories. 
The directory is determined by the state of the ITAMAS switch located in the 
common block called CQMSWT. A call to LOOK must be made before this routine 
is called to set the I.FOUND pointer. 

The calling sequence is: 

CALL PDESCR 

ITAMAS=0 Before call plant in PAGET. 

ITAMAS=1 Before call plant in XLIBD. 

Any description to be planted must be in the DESVEC vector before the call 
to the routine. DESVEC is a vector in the common block called COMDIR. A 
description of fifty characters may be planted. 

FINDING SPACE: TOP 

This routine is used to find space in the system directories. The 
directory is determined by the state of the ITAMAS switch located in the 
common block called COMSWT. A check is made to see if the upper limits of the 
directories are exceeded. The upper limits of the directories are IPAEND and 
ILDEND located in the common block called COMPNT. If the upper limit is not 
exceeded then ISPACE is set to point to the first empty location. 
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The calling statement is: 

CALL TOP 

ITAMAS=0 Before call find space in PAGET. 

ITAMAS=1 Before call find space in XLBD. 

PLANTING NAMES: PCHNAM 

This routine is similar to PLANT. Only the file name is planted. 

The calling sequence is: 

CALL PCHNAM 

PLANTING EXTENSIONS: PCHEXT 

This routine is similar to PLANT. Only the file extension is planted. 

The calling sequence is: 

CALL PCHEXT 

PLANTING PROPERTIES: PCHPRO 

This routine is similar to PLANT. Only the file property is planted. 

The calling sequence is: 

CALL PCHPRO 

CHECK FILE NAME SYNTAX: CHECKN 

This routine checks the syntax of a five or less than five character name. 
The name must be in A1 type format. A switch called ITWITC is set to indicate an 
error or correct name. LOCATT is set to the number of characters in the name. 
ITWITC and LOCATT are located in a common block called COMSWT. 

The calling sequence is: 

CALL CHECKN 

ITWITC=0 After return name error. 

ITWITC^O Correct name . 

LOCATT The number of characters in the name. 
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BLANKS should be called before the name is placed in the XPNAM buffer. 

The name must be in XPNAM before calling CHECKN. 

CONTINUATION INDICATOR: F INISH 

This routine types CONTINUE after a command has been executed. 

The calling sequence is: 

CALL FINISH 

SYSTEM DIRECTORY LISTER: LISTER 

This routine lists the system directories. The directory is determined 
by the state of the ITAMAS switch located in the common block called COMSWT. 

To list a single file name LOOK must be called to set the IFOUND pointer, 

IPASS must be set to IFOUND, LOCATT must be set to IFOUND before calling LISTER. 
To list the whole directory IPASS must be set to one and LOCATT must be set 
to zero before calling LISTER. IPASS is a variable in the common block called 
COMPAS. 

The calling sequence is: 

CALL LISTER 

ITAMAS=0 Before call list PAGET. 

ITAMAS=1 Before call list LCLIBD. 

SIXBIT DEVICE NAMES: CONDTA,CONMTA 

These routines return SIXBIT device names. 

The calling sequence for CONDTA is: 

CALL CONDTA (N, DEVICE) 

N 0< N< 6. Is one less than a Dec tape device number. 

DEVICE The returned SIXBIT device name. 

The calling sequence for CONMTA is: 
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CALL CONMTA(N, DEVICE) 

N 0< Ni 3. Is one less than a Magnetic tape device number. 



If N=3 the device is DSK. 

DEVICE The returned SIXBIT device name. 

PACKING CHARACTERS: PACK 

PACK is used to pack characters from one in a word to five in a word. 
The calling sequence is: 



CALL PACK (TEST, XSTRIN,N) 

TEST The word where the packed characters are stored. 



XSTRIN A vector of the characters to be packed. 



N 



0< N< 5. Is the number of characters to pack in TEST. 
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DISPLAY 



RAND TABLET SERVICE ROUTINE: GNF 

This is a routine that allows a user to INIT the RAND TABLET. Details 
about this routine are found in reference (2 ) . 

GENERALIZED DISPLAY: GENDIS 

This program operates the OLMGS in a general way using the RAND TABLET. 
Details about this program can be found in Part I and Part II. 

GENERALIZED RAND TABLET INPUT: GENPUT 

This program uses the OLMGS and routines GRDVAL, VALUEY, VALUEX, TPAGE, 
YESNO, RECTAN, XDATA and FINHIT to allow a generalized input graph using the 
RAND TABLET. Details about this routine can be found in Part I and Part II. 

CREATING A GENERAL GRID: GRDVAL 

This routine creates a general grid from the data limits XL, XR, YB and 
YT given as arguments . 

The calling sequence is: 



CALL GRDVAL(XL,XR,YB,YT,DC, ITEK, ISLAB,XLABEL,YLABEL) 



XL,XR 



Left and right limits of the data area. 
Bottom and top limits of the data area. 

The data area grid density. DC should never 
be a value less than 3.0; values of 8.0 to 



YB,YT 



DC 



20.0 are recommended. 



ITEK=0 

7*0 



No labels on axis . 
Label axis. 
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ISCLAB=0 

^0 



XLABEL 



YLABEL 



Standard labels. 

Scientific labels. 

A vector of fifty or less characters to be 
placed on the X-axis in A5 format. 

A vector of fifty or less characters to be 
placed on the Y-axis in A5 format. 



VALUE GENERATION FROM SCALED AREA: VALUEX, VALUE Y 

These routines use an inverse mapping on the grid scale area to obtain 
the data value of a raster coordinate. The coordinates are left in a buffer 
by the RAND TABLET service routine. Details about the RAND TABLET service 
routine can be found in reference (5) . 

The calling sequence for VALUEX is: 

CALL VALUEX (XL , XR , ICOR , VALUE) 

XL,XR Left and right limits of the data area. 

ICOR X-axis coordinate. 

VALUE The value of the scaled data area associated with 

the X-axis at this coordinate. 

The calling sequence for VALUEY is: 

CALL VALUEY (YB , YT , ICOR , VALUE) 

YB,YT Bottom and top limits of the data area. 

ICOR Y-axis coordinate. 

VALUE The value of the scaled data area associated 

with the Y-axis at this coordinate. 
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DISPLAYING CHARACTER STRINGS: TPAGE 



This routine is used to display character strings . 

The calling sequence is: 

CALL TPAGE (IX,IY,SAVF,N) 

IX, IY The raster coordinate of the point where 

the string will begin. 

SAVF A vector of characters. 

N 0 < N < 15. The number of words of characters. 

DECISION BOX: YESNO,RECTAN 

These routines generate a decision box. YESNO displays YES, NO in a box 
generated by RECTAN. 

The calling sequence for YESNO is: 

CALL YESNO (IX, I Y) 

IX, IY The coordinate of the YES. NO will be located 

at IX,IY-50. 

A box will be drawn by YESNO by calling RECTAN. 

The calling sequence for RECTAN is: 

CALL RECTAN (IX, I Y) 

IX, I Y The coordinates of the center of the box. 

GENERAL DATA DISPLAY: PDP342 

This routine uses the OLMGS to do generalized data display. Any details 
may be found in Part I . 
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The calling sequence is: 



CALL PDP342 (NPLOT , DC , NCHAP , NPTS , X , Y , I SCLAB , XLABEL , YLABEL , IDF , I ERR) 



■ NPLOT=l 


Use new graph. 


I NPLOT=2 


Use same graph. 


I DC 


Limits the density of the grid. DC should 
never have a value less than 3.0; values of 8.0 
to 20.0 are recommended. 


1 NCHAR 


0 < NCHAR < 5. The plot character selected. 
NCHAR PLOT CHARACTER 

1 

2 * 

3 X 

4 0 

5 0 


1 NPTS 


The number of points to plot. 


I X 


X- axis values . 


1 Y 


Y-axis values. 


■ ISCLAB=0 


Standard labels. 


1 ISCLAB^O 


Scientific labels. 


1 XLABEL 


A vector of fifty or less characters in A5 
format. These characters will be placed on 
the X-axis. 


I YLABEL 


A vector of fifty or less characters in ^5 format. 
These characters will be placed on the Y-axis. 


1 IDF=0 


No axis labels. 


I IDF^O 


Axis labels. 


■ IERR^l 


After return an error occurred. 
3-4 
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PARAMETER INPUT ON THE RAND TABLET: XDATA 

This routine allows a user to input data via the RAND TABLET. IRE2DI 
should be set to seme maximum data value. Details on the use of this 
routine can be found in Part I . 

The calling sequence is. 

CALL XDATA 

IRE3DI The value of the input data. 

IRE2DI and IRE3DI are located in the common block called REGPAS. 

HITTING A DECISION BOX: FINHIT 

This routine is used to flag a decision box. Details about how this 
routine works can be found in Part I. Before a call to this routine 

a call to GFNINS an entry point in GNF must be made. 

The calling statement is: 

CALL FINHIT 

IREGDI The X- coordinate of the hit. 

IRE1DI The Y-coordinate of the hit. 

IREGDI and IRE1DI are located in the common block ..called REGPAS. 

FUNCTION KEY INIT: FUNINI 

This routine is used to INIT the function keys. 

The calling sequence is: 

CALL FUNINI (ICHANN,IER) 

ICHANN 0 < ICHANN < 15. The user channel number 

IER=1 INIT error 



3-5 



O 



195 



FINDING A FUNCTION KEY DEPRESSION: KEYNUM 



This routine returns the number of the function key depressed. 

The calling sequence is: 

CALL KEYNUM (ICHANN,NUM) 

ICHANN 0 < ICHANN <15. The user channel number. 

NUM The number of the function key depressed, 

TURNING OFF FUNCTION KEY LIGHTS: LIGH TS 

This routine turns off the function key light and returns the light 
number turned off. 

The calling sequence is: 

CALL LIGHTS (ICHANN,LIGHT)_ 

ICHANN 0 < ICHANN <15. The user channel number 

LIGHT The light number turned off. 

CLOSING FUNCTION KEY CHANNEL: CLOFBX 

This routine closes the function key channel. 

The calling sequence is: 

CALL CLOFBX (ICHANN) 

ICHANN 0 < ICHANN <15. The user channel number. 

RELEASING FUNCTION KEY CHANNELS: RELFBX 

This routine releases function key channels. 

CALL RELFBX (ICHANN) 

ICHANN 0 < ICHANN <. 15. The user channel number. 

Any retails about the Function Key Monitor Service Routine may be found 
reference (7) . 
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FKLP MODE PROGRAM: FKLP 



This program assigns the OLMS commands to the function keys, 
this mode the system is operated by the function keys. 



t 
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SECTION IV 



1 INPUT-OUTPUT 

CLEAR DEC TAPE DIRECTORY: CLEAR ,CCLEAR 

This routine uses the UTPCLR program operator to clear the directory. 

The calling sequence for CLEAR is: 

CALL CLEAR (ICHANN) 

ICHANN 0 < ICHANN< 15. A user channel. 

The program CCLEAR is a CHAIN file that INITS a device and clears the 
directory by calling CLEAR. 

WRITING AND READING SYSTEM DIRECTORIES: NTAD,DIRD 

These programs are CHAIN files used to write and read the system directories 
using Dec tape and Disk. 

The use of NTAD is: 

NTAD writes a system directory on Dec tape or Disk. 

The use of DIRD is: 

DIRD reads a system directory from Dec tape or Disk. 

The device and directory is determined by the command given to INTER. 

READING AND WRITING DATA ON DEC TAPE OR DISK: GETDD , SADA 

These programs are CHAIN files used to read and write data. 

The use of GETDD is: 

GETDD reads data from Dec tape or Disk into the system buffer. 

The use of SADA is: 

SADA writes data in the system buffer on Dec tape or Disk. 

The device is determined by the command given to INTER. The system buffer is 

ERJC 
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SYSBUF located in the common block called COMBUF. 



DELETING FILES FROM DEC TAPE OR DISK: CDEDAS 

CDEDAS is a CHAIN file that deletes files from DEC tape and DISK. 

The CDEDAS program calls RENAME to delete file names. Any details about 
this program can be found in reference (1) . 



RENAMEING FILES: RENAME, CRENA 

The RENAME routine is used to rename file. 

The calling sequence for RENAME is: 

CALL RENAME (XNAME , EXT , XNAME1 , EXT1 , IER , DEVICE , 1CHANN) 



XNAME,EXT 



XNAME1,EXT1 

IER=1 

=2 

=3 

DEVICE 

ICHANN 



The name and extension of the file to be 
renamed. 

The new name and extension of the file. 

Init error 
LOOKUP error 
RENAME error 

The SIXBIT device name 

0 < ICHANN < 15. The user channel number. 



The program CRENA calls RENAME to rename files. For details see reference (1) . 



ENTERING FILE NAMES IN DEC TAPE AND DISK DIRECTORIES: ENTER 

This routine enters file names in the Dec tape and Disk directories. 
The calling sequence is: 

CALL ENTER (XNAME, EXT, IER, ICHANN) 

XNAME, EXT The file name and extension to be entered. 

IER=4 ENTER error 

ICHANN 0 < ICHANN < 15. The user channel number 
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CLOSING CHANNELS: CLOSEC 



This routine closes a channel and resets the JOBFF pointer. 
The calling sequence is: 



CALL CLOSEC (ICHANN, JOB, IER) 



JOB 



IER=1 



ICHANN 



0 < ICHANN < 15. The user channel. 
The contents of the old JOBFF pointer. 
Closing error. 



JOB and ICHANN are in the common block called INOUT. JOB is set by a call 
to BUFFER. For details see reference (1) * 

CHECKING FOR A FILE IN THE DEC TAPE AND DISK DIRECTORIES: CHECKF 

This routine does a LOOKUP on a file name. 

The calling sequence is: 



For details see reference (1) . 

WRITING AND READING MAGNETIC TAPE DIRECTORIES rMTODIR.MHDIR 

These routines are used to write and read magnetic tape directories. 

The magnetic tape directory is DIRECT. DIRECT is located in the common block 
called MTAPE. 

The calling sequences are: 



CALL CHECKF (XNAME,EXT, ICHANN, IER) 

XNAME,EXT The file name and extension to LOOKUP. 



ICHANN 



0 < ICHANN < 15. The user channel number. 



IER=2 



The file was not found. 



CALL MTODIR 



Write the tape directory. 



CALL MTIDIR 



Read the tape directory. 



t 
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POINTER ROUTINES: CDRBUF, NUMWDS 



These routines are used to manipulate pointers and plant data. 
The calling sequence for CDRBUF is: 

CALL CDRBUF (IBUF) 



IBUF 


A link variable. The contents of the right half 




of the location pointed at by IBUF is returned. 



The calling sequence for NUMWDS is 
CALL NUMWDS (IBUF, WORDS) 



IBUF 

I 


A link variable. 


WORDS 


The number to be stored. The contents of WORDS 



are stored in the right half of the location 



! 

! 


pointed at by IBUF. 


ALLOCATING BUFFERS: BUFFER 





Input or output buffers are allocated by this routine. 

The calling sequence is: 

CALL BUFFER(IND,IER, DEVICE, ICHANN, IBUF, MODE, NUMBUF, JOB) 



IND=0 


Allocate output buffers . 


IND^O 


Allocate input buffers. 


IER=1 


After return INIT error. 


DEVICE 


SIXBIT device name. 


ICHANN 


0 < ICHANN <_ 15. The user channel. 


IBUF 


A link variable that points to the buffer header 




set by the routine. 


- MODE 


The device mode. For details see reference (1) . 


r NUMBUF 


The number of buffers. 


JOB 


The contents of the old JOBFF. 
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OUTPUT: OUT I, OUT 



These routines are used to output data. 

The calling sequence for OUTI is: 

CALL OUTI (ICHANN, IER) 

ICHANN 0 <_ ICHANN <_15. The user channel number. 

IER=3 Bit transfer error. 

Before output on a channel, OUTI must be called to initialize the channel. 

The calling sequence for OUT is: 

CALL OUT (ICHANN, IER) 

ICHANN 0 < ICHANN <15. The user channel number. 

IER=4 Bit transfer error. 

OUT should be called one time for each buffer. Data may be transfered to the 
buffers using the address in the buffer header pointed at by IBUF. The routine 
called TRANSF can be used to block transfer a block of data to the output 
area. For details see reference (1). 



INPUT: INI, GET 

These routines are used to input data. 

The calling sequence for INI is: 

CALL INI (ICHANN, IER) 

ICHANN 0 < ICHANN <15. The user channel number. 

IER=6 Initialization error. 

IER=3 Bit transfer error. 

Before an input on a channel INI must be called to initialize the channel. 
The calling sequence for GET is: 



CALL GET (ICHANN, IER) 
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ICHANN 



0 < ICHANN < 15. The user channel number. 

IER=7 Input error. 

IER=5 Bit transfer error. 

IER=8 End of file was encountered. 

GET should be called one time for each buffer. Data may be transfered from 
buffers using the address in the buffer header pointed at by IBUF. The 
routine called TRANSF can be used in conjunction with GETADR to block transfer 
a whole buffer of data to the user area. For details see reference (1) • 



MAGNETIC TAPE OPERATIONS: TAPE 

This routine can do all of the magnetic tape MTAPE functions. 
The calling sequence is: 

CALL TAPE (ICHANN, FUNCTION, IER) 

ICHANN 0 < ICHANN < 15. The user channel. 

FUNCTION See reference (1) . 

IER=1 Function error. 



TESTING FILE STATUS BITS: STATUS 



This routine is used to test file status bits . 

The calling sequence is: 

CALL STATUS (ICHANN, MASS, IER) 

ICHANN 0 < ICHANN < 15. The user channel number. 

MASS An integer representing the bits to test. 

IER=1 An error occurred associated with the MASS. 

The bit information can be found in reference (1) . 



GENERAL INPUT, OUTPUT: IFILE,OFILE,IFILED,OFILED 

These routines read and write data on magnetic tape, Dec tape and Disk. 
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The calling sequence for IFILE is: 

CALL IFILE 

Before calling IFILE a pointer where the data will be stored and a pointer 
to the buffer ring must be given. ITBUF is a link variable that must point 

to the area where the data will be stored. IFBUF is a link variable that 

must point to the buffer ring. NUMBUF must be equal to the number of buffers. 
The calling sequence for OFILE is: 

CALL OFILE 

Before the call ITBUF must point to the buffer ring. IFBUF must point 

to the area where the data will be moved from during the output. NUMBUF must 

be equal to the number of buffers. 

ITBUF, IFBUF and NUMBUF are located in the common block called INOUT. 

OFILED and IFILED are defined on the same pointers as OFILE and IFILE. 

The difference is in the buffer size for DEC tapes. These routines are used 
for DEC tape. 

RELEASING CHANNELS: RELEAS 

This routine releases channels. 

The calling sequence is: 

CALL RELEAS (ICHANN, JOB) 

ICHANN 0 < ICHANN < 15. The user channel number. 

JOB The contents of the old JOBFF. 

WRITING AND READING SYSTEM DIRECTORIES ON MAGNETIC TAPE: WRTDIR,REDDIR 

These routines write and read the system directories on magnetic tape. 
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The calling sequence for WRTDIR is: 

CALL WRTDIR. 

The calling sequence for REDDIR is: 

CALL REDDIR. 

INPUTTING AND OUTPUTTING DATA ON MAGNETIC TAPE: GETF , SAVEM 

These routines input and output data on magnetic tape by calling 
IFILE and OFILE. 

The calling sequence for GETF is: 

CALL GETF. 

This call reads data. 

The calling sequence for SAVEM is: 

CALL SAVEM 

This call writes data. 



MAGNETIC TAPE CHAIN FILES: CGETF , CDIR , CNTAP , CSAVEM, CDESCR , CDELE 

These CHAIN programs are used for INPUT-OUTPUT on magnetic tape. 



CGETF 

CSAVEM 

CDIR 

CNTAP 

CDESCR 

CDELE 



Reads input data into the system buffer. 

Writes data on magnetic tape from the system buffer. 
Reads the magnetic tape directory into DIRECT. 

Writes the magnetic tape and system directories 
on magnetic tape. 

Writes file descriptions. 

Deletes file name from the magnetic tape directory.' 
This program is used to rename files on magnetic 
tape. 
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SECTION V 



MISCELLANEOUS 

SLEEPING A JOB: RESTIN 

This routine uses the SLEEP program operator to stop a job and continue 
automatically after an elapsed real time of ITIME. 

The calling sequence is: 

CALL RESTIN (ITIME) 

ITIME 0 < ITIME < 68. The number of seconds 

to sleep. 

LOGIN AND TAPE ASSIGNMENTS 
LOG 

PROJECT NUMBER { , | /^PROGRAMMER NUMBER 

Assign or Deassign Devices 

, l I l 

{AS|DEAS} 1 {<Blank>}i{<Device Name>>i 

Putting Source Deck on Disk: 

R PIP 

DSK':{<Name>}^ CDR: 

FORTRAN COMPILE 

.1 l 

COMPILE! <Blank> } j { <Name> } j 

PUTTING THE OLMS ON DISK 
R PIP 

DSK:(XB) {<Device Name> } j : { <Name> } j 
I 5-1 

ERLC 
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BUILDING A CHAIN FILE 



1 1 
LOAD {<Blank>} 1 INTOP,{<NAME>} 1 

.1 ill 

SAVE {<Blank> }j {<Device Name>} 1 {<Blank>} 1 {<Name>} 1 
RUNNING THE OLMS 

1 l l 

RUN {<Blank>}i{<Device Name>} 1 {<Blank>} 1 OLMS 

For any details see references (1) , (14) . 
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INDEX OF ROUTINES 



ROUTINE 


PAGES 


ROUTINE 


PAGES 


BLANKS 


2-3 


DELEFL 


2-1 


BUFFER 


4-4 


DIRD 


4-1 


CCLEAR 


4-1 


ENTER 


4-2 


CDEDAS 


4-2 


FINHIT 


3-5 


CDELE 


4- 8 


FINISH 


2-8 


CDESCR 


4-8 


FKLP 


3-7 


CDIR 


4-8 


FUNINI 


3-5 


CDRBUF 


4-4 


GENDIS 


3-1 


CGETF 


4-8 


GENPUT 


3-1 


CHARD 


2-3 


GET 


4-5 


CHECKF 


4-3 


GETADR 


2-3 


CHECKN 


2-7 


GETCOR 


2-1 


CLEAR 


4-1 


GETDD 


4-1 


CLOFBX 


3-6 


GETF 


4-8 


CLOSEC 


4-3 


GNF 


3-1 


CNTAP 


4-8 


GRDVAL 


3-1 


CONDTA 


2-8 


IFILE 


4-6 


CONMTA 


2-8 


IF I LED 


4-6 


CONVER 


2-4 


INI 


4-5 


CRENA 


4-2 


INTER 


2-1 


CSAVEM 


4-2 


INTER? 


2-1 
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ROUTINE 


PAGES 


ROUTINE 


PAGES 


KEYNUM 


3-6 


REDDIR 


4-7 


LIGHTS 


3-6 


RELEAS 


4-7 


LISTER 


2-8 


RELFBX 


3-6 


LISTT 


2-2 


RENAME 


4-2 


LOOK 


2-5 


RESTIN 


5-1 


MMODE 


2-3 


SADA 


4-1 


MTODIR 


4-3 


SAVEM 


4-8 


MTIDIR 


4-3 


STATUS 


4-6 


NUMWDS 


4-4 


TAPE 


4-6 


NTAD 


4-1 


TOP 


2-6 


OFILE 


4-6 


TPAGE 


3-3 


OFILED 


4-6 


TRANSF 


2-4 


OUT 


4-5 


VALUEX 


3-2 


OUTI 


4-5 


VALUEY 


3-2 


PACK 


2-9 


WRTDIR 


4-7 


PACHEXT 


2-7 


XDATA 


3-5 


PCHNAM 


2-7 


YESNO 


3-3 


PCHPRO 


2-7 






PDESCR 


2-6 






PDP342 


3-3 § 3-4 






PLANT 


2-5 






REASL 


2-9 






RECTAN 


3-3 
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